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ABSTRACT 


Anomalies  from  manufacturing  variability  or  in-service  damage  are 
manifested  as  subtle  differences  in  the  inertial  and  stiffness  matrices  of 
otherwise  identical  components.  Modal  analysis  seeks  to  model  these 
constitutive  relations  by  examining  a  structure's  dynamic  response  to  applied 
vibrations.  In  this  investigation,  characterization  of  anomalies  in  composite 
structures  is  attempted  through  multi-input,  multi-output  modal  analysis.  Two 
graphite-epoxy  specimens  manufactured  by  different  processes  are  subjected  to 
random  excitation  and  frequency  response  functions  are  obtained.  The  dynamic 
behavior  of  the  structures  is  compared  through  modal  parameter  extraction 
techniques.  Dynamic  solutions  from  a  finite  element  model  are  used  to  validate 
the  results.  Finally,  observations  are  made  concerning  the  feasibility  of 
extending  the  procedure  to  quality  assurance  and  damage  detection  applications. 
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I.  INTRODUCTION 


Traditional  methods  of  Non-Destructive  Inspection  (NDI)  for  Naval  Aviation 
applications  are  inadequate  for  identifying  damage/defects  in  composite  materials  for 
several  reasons  [Refs.  1,2],  The  widely  used  techniques  of  eddy  current,  dye 
penetrant,  and  magnetic  particle  inspection  are  applicable  only  to  metallic  structures 
and  will  only  show  structural  damage  close  to  the  surface  of  the  component. 

Ultrasonic  testing  (UT)  results  are  often  ambiguous  or  difficult  to  interpret  and 
quantify  and  require  a  reference  standard  for  calibration.  Standard  radiographic  testing 
is  difficult  to  perform  in  the  field,  presents  a  significant  safety  hazard,  is  insensitive  to 
defects  and  delaminations,  and  the  results  (in  the  form  of  film  images)  are  often 
highly  subjective  [Ref.  3].  A  technique  currently  under  development,  acoustic 
emissions  (AE)  testing,  shows  promise  but  suffers  from  sensitivity  to  background 
noise.  Clearly,  as  the  use  of  composite  materials  expands  in  programs  such  as  the  F- 
18E/F  and  V-22,  a  new  approach  to  NDI  is  needed. 

The  characterization  of  a  structure's  dynamic  behavior  may  be  achieved 
through  an  examination  of  its  response  to  an  induced  vibration,  a  process  generally 
known  as  "modal  analysis".  Modal  analysis  using  multiple  input,  multiple  output 
(MIMO)  excitation  in  the  lower  frequency  range  (below  1000  Hz)  is  relatively 
straightforward  to  perform  with  the  signal  processing  tools  now  available,  and,  given 
sufficient  resolution,  can  be  used  to  detect  differences  in  the  dynamic  response 
characteristics  of  two  seemingly  identical  structures.  The  premise  of  this  investigation 
is  that  these  differences  in  resonant  frequencies,  damping  ratios,  and  deformation 
amplitudes  can  then  be  linked  to  manufacturing  anomalies,  internal  damage,  or  other 
structural  deficiencies  in  components.  If  the  procedure  proves  viable,  the  long  term 
goal  is  to  adapt  this  technique  for  use  at  depot  and  intermediate  level  maintenance 
activities  for  both  quality  assurance  acceptance  testing  and  damage  detection. 

The  application  presented  here  sought  to  utilize  MIMO  modal  analysis  to  detect 
differences  in  the  structural  composition  of  composite  components  by  a  comparison  of 
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modal  parameters  obtained  through  vibration  tests.  Test  results  from  two  1  foot 
square  composite  plates  which  were  manufactured  by  different  processes  were 
compared  to  each  other  through  the  extraction  of  modal  parameters.  The  results  were 
validated  against  a  finite  element  model  which  can  be  updated  with  the  test  results  for 
use  as  a  baseline  in  future  measurements.  The  results  suggest  that,  though  further 
refinement  is  required,  the  procedure  is  a  viable  method  of  detecting  manufacturing 
anomalies  and  internal  damage  in  composite  materials.  Localization  of  such 
deficiencies  was  not  attempted  and  is  the  next  logical  extension  of  the  technique. 

The  damage  localization  problem  was  recently  investigated  by  Campbell 
[Ref  4],  where  frequency  domain  localization  using  the  Error  Matrix  Method  was 
applied  to  a  composite  beam  with  purposely  fabricated  damage  located  along  the 
neutral  axis.  A  single  reference  and  response  data  acquisition  system  was  used  to 
collect  the  frequency  domain  data.  In  his  investigation,  the  damage  was  consistently 
identified  but  the  localization  results  were  inconclusive.  Among  his  recommendations 
for  future  work  were  the  utilization  of  multiple  channel  data  acquisition  systems  and 
frequency  response  function  curve  fitting  techniques. 
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II.  BACKGROUND 


A.  MODAL  ANALYSIS  OVERVIEW 

Modal  analysis  refers  to  the  identification  of  the  modal  parameters  of  a 
structure  in  order  to  characterize  its  dynamic  behavior.  These  modal  parameters  are 
the  resonant  frequencies,  damping  ratios  and  mode  shapes.  The  mode  shapes  represent 
the  physical  displacement  of  the  structure  at  a  particular  resonant  frequency  (mode). 
Together,  these  quantities  make  up  a  unique  "modal  model"  of  the  structure  which 
describes  its  dynamic  behavior  and  as  such  represent  the  inertial  and  stiffness 
properties  of  the  structure.  The  modal  parameters  can  be  obtained  through  dynamic 
analysis  of  a  theoretical  (finite  element)  model  or  through  vibration  testing  of  the 
physical  structure.  Traditional  applications  of  the  modal  analysis  approach  have  been 
to  [Ref  5]: 

1.  Use  the  results  of  a  modal  test  to  verify  and  adjust  a  finite  element  model. 

2.  Produce  a  mathematical  model  of  the  structure  to  predict  the  effects  of 
structural  modifications. 

3.  Estimate  the  forces  a  structure  experiences  in  operation  through 
measurement  of  the  response  to  the  unknown  forces  and  a 
mathematical  description  of  the  transfer  function. 

The  procedure  by  which  the  modal  parameters  are  obtained  depends  upon 
whether  a  theoretical  model  or  a  vibration  test  is  used  to  derive  them.  The  theoretical 
and  experimental  approaches  to  modal  analysis  start  from  opposite  directions  to  arrive 
at  the  same  point:  identification  of  the  modal  parameters. 

B.  THEORETICAL  MODEL  APPROACH 

Modal  properties  can  be  obtained  for  a  N  Degree  of  Freedom  (DOF)  viscously 
damped  system  by  solving  the  N  equations  of  motion,  as  follows: 

[m]  +  [c]  Ig-  (t)}  +  [k]  lg(t)}  =  1 
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where 


[m]  =  mass  matrix 
[c]  =  damping  matrix 
[k]  =  stiffness  matrix 

{q(t))  =  displacement  in  generalized  coordinates 
{F(t)}  =  generalized  applied  forces 

The  system  of  equations  is  reduced  to  a  straightforward  N  by  N  eigenvalue  problem 
when  the  excitation  and  damping  are  zero  (free  vibration).  The  result  yields  the 
system's  undamped  natural  frequencies  (the  eigenvalues)  and  mode  shapes 
(eigenvectors)  for  free  vibration.  If  proportional  damping  is  assumed,  where 

tc)  -  prk]  +  Tim],  2 

then  the  mode  shapes  are  identical  to  those  derived  from  the  undamped  case  [Ref  5], 
and  the  damped  natural  frequencies  are  given  by 

to,,/  =  to/  (  1  -  C/  ),  where  =  p  to,/  2  +  y  /  2  to,  .  3 

Even  in  the  general  damped  case,  with  small  damping  ratios  the  mode  shapes  and 
frequencies  are  very  close  to  those  in  the  undamped  case. 

In  a  finite  element  model  the  total  number  of  DOF's  corresponds  to  the 
summation  of  the  DOF’s  at  each  node.  The  system  of  N  equations  can  then  be  solved 
for  the  mode  shapes  and  frequencies  using  a  number  of  algorithms  suited  for  this 
purpose  [Ref  5,  6]. 

The  response  of  the  structure  to  an  applied  excitation  is  usually  characterized  in 
the  frequency  domain  by  the  transfer  function,  known  as  the  frequency  response 
function  (  FRF  ),  from  the  applied  excitation  to  the  displacement,  velocity,  or 
acceleration  response.  The  complete  "response  model"  of  a  structure  would  consist  of 
many  different  FRF's,  one  for  each  excitation  and  response  pair.  For  a  proportionally 
damped  multiple  DOF  system,  the  FRF  can  be  written  as  [Ref.  5  ]: 
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0)^  +  i2f^0)^Ci) 


where 

A'^y  =  modal  constant  or  residue 
=  natural  frequency 
=  damping  ratio  of  r“’  mode 
N  =  total  number  of  modes  of  interest 
j  =  reference  (excitation)  coordinate 
i  =  response  coordinate 


The  coefficients  or  "residues"  are  complex  in  general  and  are  usually  expressed  in 
terms  of  amplitude  and  phase.  These  residues  represent  the  spacial  orientation  of  the 
structure  in  a  particular  mode  and  can  be  used  to  derive  the  mode  shapes.  The  FRF 
contains  information  on  all  of  the  modal  parameters,  since  the  resonant  frequencies 
occur  at  the  peaks  in  the  FRF  (w  =  w,),  and  the  relative  amplitudes  of  the  peaks 
provide  the  mode  shapes  and  damping  ratios.  Appendix  A  provides  additional 
background  on  theoretical  vibrational  analysis. 

C.  MODAL  TESTING 

In  the  modal  testing  approach,  modal  parameters  are  derived  from  examining 
the  response  of  the  system  to  a  known  excitation.  The  responses  are  typically  in  the 
form  of  accelerations,  obtained  through  accelerometers  placed  on  the  structure  at 
various  locations.  The  reference  force  measurement  is  provided  by  a  force  transducer 
connected  to  the  source  of  excitations  (shakers).  The  time  history  of  these  signals  is 
converted  into  the  frequency  domain  by  the  use  of  Fast  Fourier  Transform  (FFT) 
techniques  and  the  result  is  a  frequency  domain  transfer  function,  the  FRF,  as  H=a/f. 
As  noted  above,  the  FRF  contains  all  of  the  desired  modal  parameters;  the  challenge  is 
in  extracting  a  consistent  set  of  parameters  from  the  test  data  which  normally  consists 
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of  many  FRF's  (for  each  force  input  there  will  be  a  set  of  FRF's,  one  for  each 
accelerometer).  In  practice,  extraction  of  the  modal  parameters  from  the  set  of  FRF's 
consists  of  the  following  basic  steps: 

1 .  Obtain  FRF's  from  FFT  analyzer. 

2.  Assemble  the  set  of  FRF's  into  a  matrix  where  each  entry  in  the  matrix 
represents  the  response  at  point  i  to  an  excitation  at  point  j. 

3.  Use  a  curve-fitting  algorithm  to  obtain  the  natural  frequencies,  damping 
ratios  and  mode  shapes  for  the  functions  contained  within  the  FRF 

matrix.  An  analytic  expression  is  generated  for  the  FRF's  and  can  be 
checked  against  the  test  FRF  to  assess  the  accuracy  of  the  fit.  A  number  of 
routines  are  currently  available  for  this  purpose,  such  as  the  complex 
exponential,  direct  parameter,  and  polyreference  techniques  [Ref.  7,  8]. 

4.  Compare  the  derived  mode  shapes  through  a  correlation  method  such  as 
Modal  Assurance  Criteria  (MAC)  (Appendix  C)  to  determine  accuracy 
of  mode  shapes  and  curve  fit. 

Until  the  mid-1970's,  modal  testing  was  conducted  almost  exclusively  using  a 
single  input  force  and  often  with  only  one  response  measurement  taken  at  a  time 
[Ref.  8],  MIMO  testing  became  possible  with  the  advent  of  modem  digital  signal 
analyzers  but  the  hardware  was  relatively  expensive  and  the  parameter  estimation 
algorithms  were  not  yet  developed. 

MIMO  modal  testing  has  recently  become  more  feasible  and  inexpensive  due 
to  advances  in  piezoelectric  transducers,  digital  signal  analyzers  and  associated 
software.  MIMO  offers  several  advantages  over  single  input  excitation  methods, 
including: 

1 .  Better  distribution  of  excitation  energy,  especially  in  large  structures. 

2.  Consistent,  more  complete  FRF  database.  Measurements  are  taken 
simultaneously  vice  during  multiple  passes  for  each  spacial  location. 

3.  Ensures  all  vibrational  modes  are  adequately  excited  in  a  given  frequency 
range. 
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4.  Reduced  probability  of  introducing  systematic  errors  to  FRF's. 

These  advantages  have  made  MIMO  modal  testing  a  useful  tool  in  the  aerospace 
industry,  where  large  structures  such  as  the  Space  Shuttle  Solid  Rocket  Motors  must 
be  analyzed  [Ref  9], 

D.  TEST  ARTICLE  DESCRIPTIONS 

The  test  articles  were  two  13.750  x  14.125  x  .088  inch  plates  made  by 
Hercules  Corp  of  Magna,  Utah,  from  IM-7  graphite  epoxy  material.  Detailed  material 
property  specifications  and  photographs  of  these  plates  are  provided  in  Appendix  D. 
The  ply  stacking  sequence  and  orientation  angles  are  shown  in  Figure  1  and  were 
identical  for  each  plate;  however,  one  was  constructed  using  a  fiber  placement 
machine  (hereafter  referred  to  as  the  MLU  plate),  while  the  other  was  made  using  the 
standard  hand  lay-up  technique  (hereafter  referred  to  as  the  HLU  plate).  The  fiber 
placement  process  is  a  relatively  new  procedure  whereby  fiber  tows  are  positioned 
for  layup  by  a  special  fiber  placement  machine.  Hand  lay  up  is  the  more  conventional 
procedure  performed  by  manually  cutting  each  ply  from  B-staged  (partially  cured) 
prepreg  fabric  and  building  up  the  laminate  one  ply  at  a  time.  Like  components  made 
from  the  two  processes  are  known  to  have  different  acoustical  properties.  The  HLU 
plate  weighed  0.9520  Ibf.  and  the  MLU  plate  weighed  0.9513  Ibf. 

The  plates  are  classified  as  symmetric  angle-ply  laminates  with  generally 
orthotropic  layers  [Ref  10,  p.  165].  The  symmetric  layup  means  that  the  coupling 
stiffness  matrices  will  be  zero  and  there  will  not  be  any  coupling  between  bending  and 
extension  in  the  laminates  (see  Appendix  E).  The  extensional  and  bending  stiffness 
matrices,  however,  will  be  fully  populated  and  some  twist  coupling  stiffness  can  be 
anticipated.  This  means  that  any  curvature  from  bending  will  also  result  in  an  out  of 
plane  twist.  The  effect  of  this  twist-curvature  coupling  is  generally  to  increase  the 
strain  energy  density  and  reduce  the  vibration  frequencies.  The  effect  of  twist 
coupling  is  proportional  to  1/N,  where  N  is  the  number  of  laminate  layers  [Ref  10]. 
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CURRENT  LAM INATE 

HERCULES  IM-7/8552 
PLY  PROPERTIES 


PLY  THICKNESS  =  ,0055  INCHES 

TENSILE  STRENGTH  =  400  KSI 

SHEAR  STRENGTH  =  14,5  KSI 


TENSILE  MODULUS  =  25  MSI 

SHEAR  MODULUS  =  0,71  MSI 

FIBER  VOLUME  =  62  X 

45 

0 

-45 

90 

90 

-45 

0 

45 


Figure  1.  Test  Specimen  Ply  Layup  Sequence 
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in.  THEORETICAL  ANALYSIS 


A.  FINITE  ELEMENT  MODELING 

Finite  element  models  were  constructed  and  analyzed  using  Structural 
Dynamics  Research  Corporation  (SDRC)  IDEAS  software.  IDEAS  is  a  multipurpose 
structural  analysis  software  package  for  workstations  incorporating  design,  modeling, 
simulation,  analysis  and  testing  features.  To  verify  the  analysis  method,  finite  element 
models  were  developed  for  isotropic  aluminum  and  steel  plates  for  which  known 
natural  frequency  data  was  available.  This  enabled  the  selection  of  the  element  types 
that  would  best  predict  the  dynamic  behavior  of  the  test  specimens. 

The  model  was  built  by  first  specifying  the  geometry  and  then  using  a 
automeshing  feature  to  create  the  nodes  and  elements  at  specified  intervals.  The 
quadrilateral  parabolic  thin  shell  elements  were  selected  as  those  providing  the  most 
reasonable  results  while  being  compatible  with  the  laminate  material  feature  of  the 
software.  This  element's  bending  stiffness  is  based  on  the  Mindlin  shell  equations 
where  the  transverse  shear  strain  displacements  are  included  [Ref.  11].  Specific 
properties  of  the  parabolic  thin  shell  elements  are  available  in  Reference  1 1 .  The 
model  was  meshed  with  elements  measuring  two  inches  per  side  with  49  total 
elements  and  176  total  nodes.  Lumped  mass  elements  were  later  added  to  simulate 
the  transducers  after  their  location  was  determined.  Figure  2  shows  the  overall 
configuration. 

The  material  for  the  model  was  generated  using  the  IDEAS  Laminate  feature, 
which  allows  each  individual  ply's  material,  orientation  and  thickness  to  be  specified. 
The  plies  are  then  built  up  into  a  complete  laminate  which  was  applied  to  the  finite 
element  model.  The  material  for  our  laminate  was  the  same  for  each  ply  but  the  ply 
orientation  varied  as  noted  in  Chapter  II  above.  The  laminate  applied  to  the  finite 
element  model  was  the  same  as  that  depicted  previously  in  Figure  1. 
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Figure  2.  Finite  Element  Model  Configuration 


B.  DYNAMIC  ANALYSIS  OF  MODEL 


IDEAS  Normal  Mode  Dynamics  features  three  basic  algorithms  for  determining 
the  undamped  natural  frequencies  and  normalized  mode  shapes;  Lanzos,  Guyan 
Reduction,  and  Simultaneous  Vector  Iteration  (SVI).  Although  all  three  methods  were 
used  on  an  aluminum  beam  test  case  with  similar  results,  the  SVI  method  was  chosen 
for  the  composite  model  because  it  offers  several  advantages  noted  below  and  it  is 
compatible  with  the  IDEAS  Optimization  routine  for  model  update.  SVI  is  two  step 
inverse  power  iteration  method  that  solves  the  eigenvalue  problem 

[k]{<|)}=A.[m]{<t)}  5 

where  X  is  the  square  of  the  natural  frequencies  and  {<|)}  is  the  mode  shape  vector. 
Appendix  F  provides  details  on  the  mechanics  of  the  SVI  method. 

The  advantages  of  the  SVI  method  are  that  it  does  not  contain  the  errors 
associated  with  static  condensation,  the  iteration  can  be  restarted  for  improved 
accuracy,  and  accuracy  estimates  are  available  from  an  iteration  history  file. 

Since  the  boundary  conditions  of  our  test  approximated  that  of  a  free  structure, 
a  free  (unrestrained)  analysis  was  performed  using  "kinematic  degrees  of  freedom",  or 
KDOF.  The  use  of  KDOF's  are  required  to  account  for  the  singularities  introduced  by 
the  rigid  body  modes  of  an  unrestrained  structure.  The  rigid  body  modes  are 
eliminated  from  the  iteration  by  partitioning  the  stiffness  matrix  into  independent  and 
kinematic  DOF,  where  the  singularities  are  all  contained  in  the  KDOF  partition.  Once 
the  rigid  body  displacements  have  been  calculated  using  the  KDOF  subspace,  the 
displacements  of  the  other  DOF's  can  be  calculated  as  linear  dependents  of  the  KDOF. 
Reference  6  contains  the  details  of  this  procedure. 

In  order  for  the  rigid  body  mode  elimination  procedure  to  function,  the  number 
of  KDOF's  designated  must  be  greater  than  or  equal  to  the  number  of  rigid  body 
modes.  Six  KDOF's  were  selected  at  different  nodes  on  the  structure  such  that  if  these 
KDOF's  represented  restraints,  the  rigid  body  modes  would  be  eliminated.  The 
KDOF's  can  be  seen  as  arrows  pointing  outward  from  the  nodes  in  Figure  2.  The 
solution  was  set  up  to  estimate  the  first  20  bending  modes  and  it  converged  within 
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ten  iterations. 

C.  THEORETICAL  ANALYSIS  RESULTS 

The  analysis  identified  20  modes  in  the  0-1000  Hz  frequency  range;  the  first  16 
frequencies  are  shown  in  Table  1. 


Mode 


Freq,,  Hz 
Initial  FE 

Freq.,  Hz  FE 
w/  Trans. 

73 

63 

108 

89 

133 

112 

184 

149 

188 

157 

317 

248 

334 

251 

336 

300 

358 

321 

411 

348 

541 

403 

548 

493 

598 

1  526 

645 

551 

664 

563 

690 

579 

Table  1.  Finite  Element  Model  Solution  Frequencies 
The  mode  shapes  for  each  frequency  were  graphically  depicted  with  the  IDEAS  Post 
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Processing  feature.  Figures  3  through  5  provide  a  representative  sample.  Examination 
of  these  mode  shapes  guided  the  placement  of  the  suspension  lines  and  transducers  as 
described  in  the  next  chapter.  After  the  transducer  locations  were  determined,  the 
effect  of  these  masses  was  modeled  using  lumped  mass  elements  and  the  dynamics 
solution  was  run  again.  The  result  was  a  significant  downshift  in  the  frequencies,  as 
shown  in  Table  1.  FRF's  were  generated  for  a  unit  excitation  using  IDEAS  Model 
Response  for  the  selected  reference  and  response  locations;  a  sample  FRF  for  the  point 
inertance  (coincident  reference  and  response)  at  node  176  is  shown  in  Figure  6. 
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Figure  4.  Finite  Element  Model  Mode  Shape  (Mode  4) 
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Figure  6. 
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IV.  EXPERIMENTAL  SET-UP  AND  PROCEDURES 


A.  HARDWARE  DESCRIPTION 

The  experimental  set-up  is  depicted  in  Figures  7  through  9.  The  major 
components  are  the  test  frame,  the  shakers  and  the  data  acquisition  system,  each  of 
which  is  described  below. 

1.  Test  Frame  and  Suspension  System 

The  test  frame  was  constructed  using  1.25  inch  diameter  threaded  steel  pipe  in 
a  rectangular  shape  with  a  cross  beam  to  support  the  test  specimen.  Teflon  paste  was 
used  on  all  threaded  joints  to  provide  a  solid  connection  and  vibration  dampening. 
Floor  flanges  attached  to  each  leg  supported  the  frame  on  a  wooden  workbench. 
Rubber  matting  was  placed  under  each  leg  to  provide  some  vibration  isolation  from 
the  shakers.  Provisions  were  made  to  attach  guy  wires  to  the  upper  comers  of  the 
structure  in  the  event  that  further  rigidity  was  necessary.  Two  insulated  clamps  and  a 
series  of  steel  tumbuckles  were  connected  to  elastic  bands  and  monofilament  line 
to  suspend  the  test  specimen  from  the  cross  beam.  This  system  provided  a  nearly 
"free"  test  condition  in  that  the  rigid  body  modes  of  the  test  article  were  at 
frequencies  below  20  %  of  the  first  bending  mode  [Ref.  5].  To  ensure  minimum 
interference  from  the  rigid  body  modes  on  the  first  bending  mode,  the  suspension 
lines  were  positioned  as  close  to  the  nodes  of  the  first  bending  mode  as  possible,  as 
determined  by  the  FE  analysis.  This  position  corresponded  to  approximately  3  inches 
from  the  edges  of  the  plate.  As  a  further  measure  in  achieving  the  free  condition,  the 
test  article  was  suspended  in  a  vertical  orientation  so  that  the  suspension  lines  were 
peipendicular  to  the  primary  direction  of  vibration. 
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Component  Description 

1.  Hewlett-Packard  715/75  Model  A2896A  Workstation 

2.  HP  A2094A  19"  Color  Monitor 

3.  HP  C1429B  Keyboard 

4.  10-Base  2  Ethernet  Transceiver 

5.  Zonic  System  7000  Analog  Signal  Chassis  Model  7014 

6.  Zonic  System  7000  Digital  Signal  Chassis  Model  7015 

7.  Techron  Model  7521  Power  Amplifier  (2) 

8.  Ling  Electronics  Model  207  Shaker  (2) 

9.  PCB  ICP  Model  208B  Force  Transducer  (2) 

10.  PCB  ICP  Model  336C34  Acceleration  Transducer  (14) 

11.  PCB  Model  070C29  Patch  Panel 


Figure  7.  Experimental  Set-Up 


2.  Shakers  and  Power  Ampliflers 

Excitation  force  was  provided  using  two  identical  Ling  Electronics  Model  207 
electrodynamic  shakers  with  an  rms  force  output  of  5  Ibf  in  the  random  mode,  a  90  g 
peak  acceleration,  and  a  frequency  range  of  DC- 10,000  Hz.  The  shakers'  resonant 
frequency  was  approximately  7.5-8  kHz,  far  outside  the  range  of  our  application.  The 
shakers  were  mounted  in  a  trunnion  frame  that  allowed  the  shaker  armature  mounting 
surface  to  be  swiveled  90  degrees  to  the  horizontal  position.  Force  was  applied 
using  "stingers"  made  from  3  inch  threaded  steel  rods  which  were  inserted  into  the 
internal  threads  of  the  armature  mounting  surface.  The  free  end  of  the  stingers  were 
screwed  directly  into  the  threaded  receptacles  of  the  force  transducers. 

The  shakers  were  each  powered  by  separate  Techron  Model  7521  amplifiers 
capable  of  delivering  a  maximum  voltage  gain  of  26.3  +/-  0.3  dB  at  150  watts  (2 
ohm  load)  and  with  a  frequency  response  of  DC  to  35  kHz.  The  amplifiers  received 
the  excitation  signals  from  the  output  module  of  Zonic  System  7000  and  applied  the 
amplified  signals  directly  to  the  shakers. 

3.  Transducers 

Force  and  acceleration  measurements  were  acquired  through  the  use  of  two 
PCB  ICP  Model  208B  quartz  force  transducers  and  14  PCB  ICP  Model  336C34  quartz 
acceleration  transducers.  The  force  transducers  had  a  working  range  of  0-10  Ibf  and 
.01-14,000  Hz  and  a  sensitivity  of  500  mV/lbf  The  accelerometers  had  working 
ranges  of  0-500  g's  and  2-5000  Hz  and  a  sensitivity  of  10  mV/g.  Calibration  of  these 
sensors  had  been  accomplished  by  PCB  within  six  months  of  the  experiment  and 
calibration  data  was  provided  for  each  transducer.  The  ICP  (Integrated  Circuit 
Piezoelectric)  sensors  use  the  piezoelectric  effect  of  a  quartz  crystal  to  generate  a 
charge  signal  proportional  to  the  applied  force  or  acceleration.  The  PCB  sensors  have 
built-in  signal  conditioning  electronics  to  convert  the  high  impedance  piezoelectric 
signal  to  a  low  impedance  voltage  signal  that  can  be  effectively  transmitted  to  the 
signal  analyser.  The  transducers  require  a  constant  20  VDC  from  the  signal  analyser 
to  operate. 
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The  accelerometers  were  attached  to  the  test  article  using  Petrowax,  which 
provides  adequate  adhesion  with  a  minimum  of  output  distortion  in  the  frequency 
range  of  interest  (less  than  2  %  below  1000  Hz).  The  use  of  Petrowax  allowed  easy 
removal  and  relocation  of  the  accelerometers.  The  accelerometers  were  connected  to  a 
patch  panel  via  lightweight  cable  with  modular  connections  at  both  ends.  The  patch 
panel  was  connected  to  the  Zonic  System  7000  via  a  16  channel  coaxial  cable  with 
separate  BNC  connectors.  The  force  transducers  were  mounted  directly  to  the  test 
specimens  using  double  sided  adhesive  tape.  The  force  transducers  were  connected  to 
the  patch  panel  via  general  purpose  coaxial  cable. 

4.  Signal  Analyser 

A  Zonic  System  7000  signal  processor  provided  multi-channel  data  acquisition, 
real  time  Fast  Fourier  Transform  signal  analysis,  and  source  signal  generation.  The 
System  7000  consisted  of  an  analog  chassis  and  a  digital  chassis,  with  the  following 
functions. 

a.  Analog  Chassis 

The  analog  chassis  contained  the  Analog  Signal  Conditioning  Input 
(ASC-I)  module  and  the  Analog  Signal  Output  (ASC-0)  module.  The  ASC-I 
performed  the  data  acquisition  including  anti-alias  filtering,  overload  detection, 
amplitude  range  selection  and  digitizing  of  the  analog  input  signals.  It  processed  16 
input  channels  up  to  5  KHz,  with  separate  A  to  D  conversion  for  each  channel  in  a 
parallel  arrangement.  Filtering  was  accomplished  automatically  as  a  function  of  the 
sampling  rate,  and  was  identical  for  all  16  channels.  The  digitized,  differentially 
balanced  signal  was  passed  to  the  digital  chassis  via  a  25  pin  RS-422  connector. 

The  ASC-O  received  digital  source  signals  from  the  digital  chassis  and 
converted  them  to  analog  signals  which  were  output  through  BNC  connectors.  The 
model  used  in  the  experiment  had  a  maximum  of  two  output  channels  available.  The 
analog  output  signals  were  sent  through  the  amplifiers  and  on  to  the  shakers  to 
produce  the  excitation.  The  system  was  capable  of  generating  sine,  random,  pseudo 
random,  and  burst  random  output  signals. 
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b.  Digital  Chassis 

The  digital  chassis  contained  the  Digital  Signal  Conditioning  Input 
(DSC-I)  module,  the  Digital  Signal  Conditioning  Output  (DSC-O)  module  and  the 
Digital  Signal  Processing  (DSP)  module.  The  DSC-I  provided  digital  filtering, 
buffering,  gain/offset  management,  zoom  processing,  trigger  control  and  ASC-I  control 
functions  via  integer  array  processors.  It  provided  data  in  a  parallel  digital  format  to 
the  DSP  for  processing. 

The  DSC-O  generated  specified  waveforms  in  a  digital  format  and  fed 
them  to  the  ASC-0  where  they  were  converted  to  analog  signals.  It  also  provided 
digital  filtering,  bursting  and  ramping  of  the  output  signals  and  interfaced  with  the 
DSP  for  control  of  the  ASC-0. 

The  DSP  provided  windowing,  signal  averaging,  trigger  control  and 
performed  variable  block  size  FFT  for  each  input  channel.  It  interfaced  with  the  host 
computer  via  an  Ethernet  Nodem  and  controlled  the  functions  of  the  DSC's. 

5.  Host  Computer 

The  data  acquisition  system  was  controlled  by  a  Hewlett-Packard  715/75  Model 
A2896A  Workstation  configured  with  two  1  GB  hard  drives  and  a  19  inch  color 
monitor.  The  computer  communicated  with  the  System  7000  signal  analyser  via  a  10- 
base-2  ethemet  cable  and  transceiver.  HP-UX,  the  Hewlett-Packard  version  of  Unix, 
and  HPVUE  (XWindows)  were  used  as  the  operating  system  on  which  Structural 
Dynamics  Research  Corporation  (SDRC)  IDEAS  was  loaded.  IDEAS  is  a  structural 
analysis  software  application  capable  of  performing  design,  drafting,  finite  element 
modeling,  dynamic  analysis,  and  testing.  It  was  used  to  control  the  experiment, 
acquire  the  data,  and  conduct  post  test  analysis.  Test  data  and  finite  element  analysis 
results  were  stored  on  one  of  the  two  internal  hard  drives.  A  DDS  SCSI-2  tape  drive 
was  used  to  back  up  the  data  collection. 
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B.  PRE-TEST  ANALYSIS 


A  number  of  trial  tests  were  conducted  using  isotropic  aluminum  and  plexiglass 
plates  to  determine  the  optimum  choice  of  data  acquisition  parameters,  excitation 
signals,  and  test  article  suspension  configuration. 

1.  Excitation  Signal 

Various  forms  of  excitation  signals  were  available  with  the  hardware  and 
software  configuration,  including  sine,  pure  random,  burst  random  and  periodic 
random  (swept  sine  is  not  available  for  multiple  input  configurations  in  IDEAS). 

Based  on  Reference  12,  which  provides  a  survey  of  these  types  of  signals,  and  some 
experimentation  with  the  different  signals,  burst  random  excitation  was  chosen.  It 
provides  excitation  in  all  frequency  ranges  of  interest  and  is  much  less  time  consuming 
than  a  manual  sine  sweep  of  the  frequencies.  There  was  very  little  difference  in  the 
FRF's  generated  from  a  single  input  swept  sine  and  the  random  excitation,  provided 
measurement  frame  averaging  was  utilized.  Figure  10  is  a  comparison  between  FRF's 
produced  by  swept  sine  and  random  excitation  for  a  plexiglass  beam,  where  the  swept 
sine  resolution  was  10  Hz.  Burst  random  excitation  was  utilized  to  combat  leakage 
problems  since  the  signal  is  brought  to  zero  well  before  the  end  of  the  measurement 
frame. 

As  discussed  in  references  5  and  12,  it  is  important  to  have  uncorrelated 
excitation  signals  when  using  multiple  inputs.  This  is  due  to  the  procedure  by  which 
the  FRF's  are  derived  from  random  excitation.  Appendix  H  gives  the  mathematical 
development,  but  suffice  it  to  say  here  that  the  FRF's  are  generated  from  the  transfer 
functions  between  the  input  and  response  spectral  densities,  cross  spectrums,  or  both. 

If  the  off  diagonal  terms  of  the  input  cross  spectrum  are  not  zero  (indicating  correlated 
inputs),  inaccuracies  in  the  FRF  calculations  can  result.  The  Zonic  signal  analyser 
provided  two  uncorrelated  random  source  signals  that  met  this  criteria. 

2.  Data  Acquisition  Parameters 

The  IDEAS  Testing  and  Zonic  System  7000  configuration  offers  a  wide  variety 
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of  signal  acquisition  and  processing  options,  including  anti-aliassing,  windowing, 
frame  averaging  and  FFT  analysis  methods.  After  some  experimentation  with  various 
combinations  of  the  available  acquisition  parameters,  the  following  settings  were  found 
to  produce  the  best  combination  of  high  signal  to  noise  ratio  and  reproducibility. 
Appendix  I  provides  the  specific  acquisition  settings  for  the  test  runs  used  in  the 
analysis. 

a.  Windowing 

A  Hanning  Broad  window  provided  low  leakage  and  good  frequency 

definition. 

b.  Averaging 

Frame  averaging  with  automatic  rejection  of  frames  with  out  of  range 
signals  was  utilized;  20  frames  were  averaged  for  each  acquisition. 

c.  Spectral  Lines  /  Sampling  Rate 

801  spectral  lines  were  used  over  the  1000  Hz  frequency  range  of 
interest,  corresponding  to  a  frequency  resolution  of  1.25  Hz.  The  sampling  rate  was 
fixed  by  the  number  of  spectral  lines  and  the  Niquist  criteria  to  2560  samples  per 
second  with  a  frame  length  of  0.8  seconds.  Total  test  time  with  20  frames  was  then 
16  seconds. 

3.  Sensor  Location 

The  finite  element  model  analysis  provided  the  means  to  determine  the  best 
positions  for  the  sensors.  Essentially,  the  exciter  locations  were  chosen  such  that  they 
were  not  coincident  with  any  of  the  nodes  predicted  by  the  FE  mode  shapes. 
Accelerometer  locations  were  arranged  to  give  a  good  spacial  representation  of  the 
structure  by  spreading  them  out  over  the  surface  without  regard  to  symmetry. 

Locations  of  both  force  and  acceleration  transducers  are  shown  in  Figure  1 1 . 

4.  Frequency  Range 

The  frequency  range  on  which  to  focus  our  analysis  was  also  determined  by 
the  finite  element  model.  The  desired  characteristics  were  well  defined  modes  that 
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Figure  1 1 .  Transducer  Locations  and  Channel  Assignment 
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were  not  unduly  influenced  by  residual  stiffness  or  inertia.  Rigid  body  modes  were 
not  of  interest  so  that  a  low  frequency  cut-off  was  established  at  20  Hz.  By 
examining  the  FRF's  generated  from  the  theoretical  analysis,  it  was  decided  that  the 
first  16  modes  met  the  above  criteria  and  the  upper  frequency  cut-off  was  therefore  set 
at  1000  Hz. 

C.  TEST  PROCEDURES 

1.  Calibration 

Prior  to  conducting  the  experiment  the  system  was  run  through  a  calibration 
sequence  on  two  levels.  First,  the  Zonic  System  7000  was  calibrated  internally  using 
a  standard  routine  from  the  IDEAS  Test  interface  that  uses  an  internal  calibration 
signal  to  set  gain  control  and  phase  offset.  Second,  the  individual  transducers  were 
calibrated  by  setting  scale  factors  in  a  mass  calibration  technique  described  in 
Appendix  J. 

2.  Data  Collection 

After  calibration,  the  transducers  were  connected  to  the  test  article  as  described 
in  Section  A  above.  The  experiment  was  controlled  from  the  IDEAS  Test  interface, 
which  allowed  for  various  options  on  how  the  data  was  collected  and  stored.  Of  the  16 
available  input  channels,  two  were  designated  as  reference  (force  input)  channels  and 
the  remaining  14  were  designated  as  response  (acceleration)  channels.  The  IDEAS 
input  channel  configuration  scheme  allowed  the  individual  transducer  sensitivities  to 
be  set  for  each  channel.  Data  acquisition  parameters  were  set  as  discussed  in  Section 
B  above  and  are  listed  in  Appendix  I.  The  source  signals  were  started  and  data  was 
acquired,  with  the  actual  data  collection  lasting  approximately  16  seconds.  Real  time 
signals  were  displayed  for  all  16  channels  during  the  test.  Time  histories  of  the  data 
were  not  stored;  rather,  the  FRF's  for  each  reference-response  pair  (28  total)  were 
automatically  computed  and  available  for  display  and  storage. 
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3.  Initial  Reproducibility  Assessment 

Multiple  runs  for  each  test  specimen  were  conducted,  stored  and  compared  to 
ensure  reproducibility  (  1 1  tests  were  stored  for  the  HLU,  9  for  the  MLU  ).  Figure 
12  shows  an  overlay  of  two  FRF's  for  the  HLU  from  different  tests  and  shows  the 
degree  of  reproducibility  obtained.  Figure  13  is  a  similar  overlay  from  two  different 
MLU  tests.  In  this  case,  all  transducers  had  been  removed  and  the  plate  had  been 
disconnected  from  the  test  fixture  between  tests.  (Both  Figures  12  and  13  are  point 
inertances  for  node  52,  where  "node"  refers  to  the  FE  model  node.  (See  Figure  1 1  for 
location.)  The  frequencies  at  all  modes  appear  to  correlate  very  well;  however, 
differences  in  amplitude  were  noted  for  some  modes  as  shown  in  both  Figures  12  and 
13.  These  differences  in  amplitude  represent  the  baseline  error  for  the  experiment.  A 
more  quantitative  error  assessment  is  made  in  the  following  section. 

Tests  were  run  using  two  different  methods  for  computing  the  FRF's  (H,  and 
H^,  as  discussed  in  Appendix  H)  and  compared,  with  little  discemable  difference  noted 
in  amplitude  and  sharpness  of  peaks,  as  shown  in  Figure  14.  In  general,  the  H,  FRF's 
appeared  to  be  of  better  quality  (less  noise)  and  were  used  for  extraction  of  modal 
parameters. 

D.  RESULTS  AND  ANALYSIS 

1.  General  Comparison  Between  Specimens 

With  sets  of  28  FRF's  for  each  test  collected  and  stored,  an  analysis  of  the  data 
commenced.  Figure  15  shows  an  overlay  of  the  point  inertance  FRF's  for  node  176 
for  each  plate.  Figure  16  is  a  similar  plot  for  node  52.  Because  of  the  density  of 
modes  and  the  fact  that  the  higher  frequency  modes  were  less  well  defined,  the 
analysis  was  confined  to  the  20-600  Hz  range,  in  which  approximately  16  well  defined 
modes  can  be  discerned.  The  overlays  in  Figures  15  and  16  indicate  an  appreciable 
difference  between  the  two  FRF's  in  amplitude,  phase,  and  frequency.  In  order  to 
quantify  and  substantiate  these  differences  over  all  reference  locations,  the  analysis 
proceeded  to  extract  the  modal  parameters  of  each  plate. 
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Figure  13.  MLU  FRF  Reproducibilit) 


33 


c  t  1  o  n 

- ^ j ^ 


4  :  S2Z*  52Z**  1 

E  TEST  9  /  20  MAR  95 


Node  52 


Phase 


Acce I ero t i on 


20  200  400  599 


Frequency  <Hz> 

2  ;176Z*  1762-  1  2:1 76 2+  1762-1 

LU  PLPTE  TEST  21  /  04  MPR  95  MLU  PLATE  TEST  2  /  05  MAR  95 


Figure  15.  HLU/MLU  FRP  Comparison  ,  Node  176 
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2.  Extraction  of  Modal  Parameters 

The  extraction  of  the  modal  parameters  was  accomplished  entirely  within 
IDEAS  Modal  Analysis  and  involved  a  number  of  steps  including  using  a  Multivariate 
Mode  Indicator  Function  (MMIF);  building  the  FRF  matrix;  estimating  the  natural 
frequencies,  damping  ratios,  and  modal  constants  (residues);  generating  the  analytic 
functions  that  model  the  FRF’s;  and,  finally,  calculating  a  set  of  mode  shapes.  The 
Polyreference  curve  fitting  technique  was  utilized  to  estimate  the  natural  frequencies 
and  damping  ratios  and  the  Frequency  Polyreference  Residue  technique  was  used  to 
calculate  the  residues  and  mode  shapes.  The  Polyreference  techniques  obtain 
global  least  squares  estimates  of  the  modal  parameters  by  processing  all  the  FRF's 
simultaneously.  The  Frequency  Polyreference  Residue  technique  has  the  advantage  of 
operating  in  the  frequency  domain  so  that  residual  inertial  and  compliance  terms  can 
be  added  to  account  for  the  effect  of  modes  outside  the  frequency  range  of  interest.  A 
more  detailed  explanation  of  the  procedure  is  presented  in  Appendix  B. 

The  end  result  of  the  extraction  procedure  is  a  parameter  table  listing  the 
natural  frequencies,  damping  ratios,  and  residue  amplitudes  and  phases  for  each 
identified  mode.  An  analytic  expression  for  the  FRF  is  generated  from  these 
parameters  to  check  the  fit  with  the  measured  FRF.  Figures  17  and  18  show  the 
analytically  derived  functions  and  the  measured  counterparts  at  node  52  for  the  HLU 
and  MLU,  respectively. 

Included  in  the  parameter  table  is  a  Modal  Confidence  Factor  (MCF),  which  is 
a  numerical  value  from  0  to  1  that  can  be  used  to  discriminate  between  computational 
(false)  modes  and  physical  (true)  modes.  Computational  modes  are  generated  by  the 
curve  fitting  algorithm  to  account  for  unwanted  effects  of  noise,  leakage,  and 
nonlinearities  and  do  not  represent  physical  modes.  An  MCF  close  to  1  indicates  a 
physical  mode,  while  an  MCF  less  than  one  indicates  a  computational  mode.  This 
feature  allows  easy  filtering  of  computational  modes  by  setting  a  MCF  threshold  value 
for  the  parameter  table.  Setting  an  MCF  threshold  of  0.95  resulted  in  16  modes  in  the 
0-600  Hz  region  for  both  HLU  and  MLU  tests.  The  resulting  parameter  tables  for  the 
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Figure  18.  MLU  Analytic  FRF  (Dashed)  and  Measured  FRF  (Solid),  Node  52 
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HLU  and  MLU  tests  are  listed  in  Appendix  J.  Table  2  shows  a  comparison  of  the 
HLU/MLU  frequencies  where  it  is  evident  that  14  out  of  16  MLU  frequencies  are 
higher.  (The  number  listed  after  "HLU"  and  "MLU"  in  the  table  indicates  the  number 
of  the  test  run.)  This  frequency  shift  is  portrayed  graphically  in  Figure  19.  The  FE 
model  frequencies  are  also  listed  in  Table  2  for  comparison. 

A  similar  comparison  for  the  damping  ratios  (in  percent  of  critical  damping)  is 
shown  in  Table  3  and  Figure  20,  where  all  of  the  MLU  damping  ratios  are  higher. 

Note  also  that  many  of  the  damping  ratios  are  quite  high  (greater  than  one  percent), 
particularly  for  the  second  mode.  This  could  be  due  to  insufficient  spectral  resolution 
(true  maxima  not  identified),  non-proportional  damping,  or  damping  introduced  by  the 
support  conditions. 

3.  Mode  Shape  Correlation 

A  Modal  Assurance  Criteria  (MAC)  matrix  was  calculated  in  order  to  correlate 
the  mode  shapes.  The  MAC  method  uses  a  least  squares  technique  to  correlate  two 
sets  of  mode  shapes,  either  from  two  different  reference  coordinates  of  the  same 
structure,  two  different  tests,  or  from  a  test  and  a  finite  element  model  solution.  A 
MAC  value  of  one  indicates  a  perfect  correlation,  whereas  a  MAC  value  near  zero 
indicates  linear  independence  between  mode  shapes.  When  assembled  into  a  matrix, 
the  rows  are  the  modes  (or  DOF's)  for  one  set  of  mode  shapes  and  the  columns  are 
the  modes  for  the  other  set  of  mode  shapes.  The  MAC  matrix  is  typically  presented 
graphically  with  the  rows  and  columns  of  the  matrix  making  up  the  horizontal  floor 
and  the  MAC  value  making  up  the  vertical  bars.  Appendix  C  has  more  detailed 
information  on  calculating  the  MAC. 

The  first  correlations  were  performed  between  mode  shapes  calculated  from 
different  reference  locations  on  the  same  plate  to  check  the  adequacy  of  the  data  and 
curve  fit.  Figures  21  and  22  show  the  MAC  matrices  for  the  HLU  and  MLU  plates, 
respectively,  where  a  high  correlation  between  mode  shapes  is  evident  from  the  MAC 
values  close  to  1  along  the  diagonal.  The  non-zero  values  for  off-diagonal  MAC 
terms  are  probably  due  to  some  noise  in  the  data.  The  MAC  matrix  for  the  HLU  / 
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1  MODE 

HLU21  FREQ 

MLU2  FREO 

DIFFERENCE 

%  DIFF 

FEFREO 

64.084 

67.955 

3.871 

6.04 

63.12 

82.878 

86.887 

4.009 

4.84 

88.5 

110.628 

114.458 

3.83 

3.46 

112.27 

133.179 

138.645 

5.466 

4.10 

148.74 

159.559 

163.659 

4.1 

2.57 

157.16 

221.542 

221.525 

-0.017 

0.01 

247.84 

238.048 

242,358 

4.31 

1.81 

251.21 

8 

288.421 

291.394 

2.973 

1.03 

300.4 

■■ 

306.315 

304.421 

-1.894 

0.62 

321.46 

343.383 

346.639 

3.256 

0.95 

347.88 

390.374 

392.949 

2.575 

0.66 

403.1 

465.919 

468.634 

2.715 

0.58 

— 

■  ■ 

494.100 

500.195 

6.095 

1.23 

493.75 

505.978 

521.241 

15.263 

3.02 

526.34 

n 

553.979 

558.03 

4.051 

0.73 

551.64 

■Bi 

562.869 

565.273 

2.404 

0.43 

563.87 

Table  2.  HLU/MLU  Frequency  Comparison 


Figure  19.  HLU/MLU  Frequency  Difference  Chart 
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HLU21  DAMPING 

MLU2  DAMPING 

DIFFERENCE 

0.908 

3.081 

2.173 

239.317 

4.054 

4.407 

0.353 

8.707 

0.852 

1.621 

0.769 

90.258 

0.514 

1.591 

1.077 

209.533 

0.215 

1.318 

1.103 

513.023 

0.587 

0.746 

0.159 

27.087 

0.189 

0.902 

0.713 

377.249 

0.412 

0.810 

0.398 

96.602 

0.399 

0.925 

0.526 

131.830 

0.227 

0.975 

0.748 

329.515 

0.781 

2.480 

1.699 

217.542 

0.579 

0.905 

0.326 

56.304 

0.776 

0.906 

0.130 

16.753 

0.510 

1.098 

0.588 

115.294 

0.299 

0.621 

0.322 

107,692 

0.436 

0.920 

0.484 

111.009 

Tables.  HI 

.U/MLU  Dampinj 

%  Ratio  Comparis( 

on 

Figure  20.  HLU/MLU  Damping  Ratio  Difference  Chart 
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Figure  22.  MLU  Modal  Assurance  Criteria  Matrix 
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MLU  correlation  is  shown  in  Figures  23  and  24  for  the  real  and  complex  mode 
shapes,  respectively.  It  is  evident  from  these  figures  that  poor  correlation  exists  for  a 
number  of  modes,  as  seen  by  the  low  MAC  values  along  the  diagonal.  Appendix  J 
lists  the  diagonal  MAC  values  for  the  complex  shapes. 

As  an  additional  comparison  between  the  mode  shapes  of  the  plates,  the 
IDEAS  Post  Processing  feature  was  utilized  to  graphically  portray  the  physical 
deflections  of  the  two  plates.  Figures  25  through  30  show  a  comparison  between  the 
two  plates  and  the  FE  model  at  modes  1,  13  and  16,  where  a  qualitative  3-D  view  is 
shown  first  and  a  side  view  for  a  more  quantitative  assessment  is  shown  second. 

(Note  that  only  the  nodes  coincident  with  an  accelerometer  position  are  shown  as 
deflected  for  the  experimental  data.)  At  the  first  mode  (Figures  25  and  26),  a  good 
correlation  exists  (MAC  =  .98)  and  the  deflections  closely  agree.  For  mode  13 
(Figures  27  and  28,  MAC=.723),  a  significant  difference  can  be  seen  in  the 
deflections  at  the  upper  left  and  lower  right  portions  of  the  plates  (MAC=.056). 
Finally,  Figures  29  and  30  show  similar  results  for  mode  16,  where  differences  are 
seen  in  the  upper  right  and  left  comers  of  the  plates. 

Before  drawing  any  conclusions  about  the  meaning  of  the  apparent  difference 
in  mode  shapes,  it  should  be  noted  that  there  are  several  possible  explanations  for  the 
poor  correlation  [Ref  7,  p.  677]: 

1 .  The  tested  system  is  nonstationary  due  to  changes  in  mass,  stiffness,  or 

damping  properties  during  the  test. 

2.  The  system  is  non-linear. 

3.  The  reference  mode  shape  contains  noise. 

4.  The  method  used  for  parameter  extraction  is  invalid  for  the  measured  data. 

5.  The  mode  shapes  are  linearly  independent. 

If  the  first  four  of  these  reasons  can  be  eliminated,  the  differences  in  mode  shapes  can 
be  considered  to  represent  a  real  difference  in  the  dynamic  behavior  of  the  two  plates. 
The  first  possibility  can  be  discarded  by  observing  that  the  test  environment 
(temperature,  pressure  suspension  method,  etc.)  did  not  change  appreciably  between 
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Figure  23.  HLU/MLU  Modal  Assurance  Criteria  Matrix,  Real  Modes 
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Figure  24.  HLU/MLU  Modal  Assurance  Criteria  Matrix,  Complex  Modes 
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Figure  26.  Mode  Shape  Comparison,  Side  View,  Mode  1 
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Figure  27.  Mode  Shape 


parison.  Qualitative  View,  Mode  13 
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Figure  28.  Mode  Shape  Comparison,  Side  View,  Mode  13 
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Figure  29.  Mode  Shape  Comparison,  Qualitative  View,  Mode  16 


Figure  30.  Mode  Shape  Comparison,  Side  View,  Mode  16 


tests  so  that  no  changes  in  material  properties  would  have  occured.  Secondly,  the 
forces  applied  to  the  structure  were  not  high  enough  to  cause  non-linear  behavior  in 
the  test  specimens.  Examination  of  the  FRF's  does  not  indicate  an  inordinate  amount 
of  noise  and  the  fact  that  the  tests  were  reproducible  supports  this  observation. 
Concerning  the  fourth  reason  above,  the  natural  frequencies  derived  by  the 
polyreference  technique  correlate  fairly  well  with  the  Finite  Element  solution,  as 
shown  earlier  in  Table  2,  substantiating  the  extraction  results. 

4.  Error  Estimation 

The  minimum  frequency  error  is  defined  by  the  spectral  resolution,  which  in 
the  case  of  these  measurements  was  1.25  Hz.  To  better  assess  the  true  error  and  the 
degree  of  reproducibility,  modal  parameters  were  extracted  from  several  different 
MLU  test  runs  and  compared  against  the  HLU  parameters  obtained  previously.  Tables 
4  and  5  and  Figures  30  and  3 1  show  the  frequency  comparisons  which  indicate  a 
similar  shift  to  that  observed  in  Figure  19,  though  it  is  somewhat  less  pronounced. 
Tables  6  and  7  and  Figures  32  and  33  show  a  comparison  between  frequencies  derived 
from  different  MLU  test  runs.  Figure  32  is  a  comparison  of  two  consecutive  runs 
conducted  one  day  apart  where  the  transducers  remained  in  place.  Figure  33  is  a 
comparison  of  two  MLU  tests  (runs  2  and  9)  where  all  of  the  transducers  had  been 
removed  and  the  plate  unhooked  from  the  support  frame  between  runs.  The  reduction 
in  reproducibility  when  the  transducers  are  removed  and  replaced  is  obvious  from 
these  figures. 

A  similar  comparison  of  the  extracted  damping  ratios  was  carried  out  and  is 
illustrated  in  Tables  8  through  11  and  Figures  35  through  38.  The  bias  towards  higher 
damping  ratios  for  the  MLU  seen  in  Figure  20  is  evident  in  Figures  35  and  36; 
however,  the  comparison  between  MLU  tests  in  Figures  37  and  38  shows  that 
significant  variation  in  damping  ratios  exists  between  test  runs  on  the  same  plate. 

The  fact  that  a  noticeable  bias  exists  between  MLU  test  runs  2  and  9  (Figure  38), 
where  transducers  were  removed  and  replaced,  and  not  for  test  runs  8  and  9  (Figure 
37),  where  transducers  remained  in  place,  indicates  that  the  damping  bias  may  have 
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Figure  32.  HLU  Run  21/MLU  Run  9  Frequency  Difference  Chart 
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MODE 

MLU9  FREO 

MLU8  FREO 

DIFFERENCE 

%  DIFF 

1 

66.127 

66.080 

-0.047 

0.07 

85.742 

82.250 

-3.492 

4.07 

112.608 

112.008 

-0.6 

0.53 

140.719 

138.099 

-2.62 

1.86 

160.617 

161.354 

0.737 

0.46 

223.014 

222.916 

-0.098 

0,04 

243.535 

243.240 

-0.295 

0.12 

8 

291.771 

291.783 

0.012 

0.00 

■■ 

304.634 

304.379 

-0.255 

0.08 

349.821 

349.264 

-0.557 

0.16 

■H 

398.932 

396.721 

-2.211 

0.55 

468.320 

469.701 

1.381 

0.29 

503.548 

502.884 

-0.664 

0.13 

522.286 

520.150 

-2.136 

0.41 

560.423 

561.094 

0.671 

0.12 

567.995 

567.520 

-0.475 

0.08 

Table  6.  MLU  Runs  8  and  9  Frequency  Comparison  (Reproducibility) 


Figure  33.  MLU  Runs  8  and  9  Frequency  Difference  Chart  (Reproducibility) 
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MODE 

MLU9FRE0 

MLU2  FREQ 

DIFFERENCE 

%  DIFF 

1 

67.955 

1.828 

2.76 

86.887 

1.145 

1.34 

112.608 

114.458 

1.85 

1.64 

140.719 

138.645 

-2.074 

1.47 

160.617 

163.659 

3.042 

1.89 

223.014 

221.525 

-1.489 

0.67 

243.535 

242.358 

-1.177 

0.48 

8 

291.771 

291.394 

-0.377 

0.13 

MM 

304.634 

304.421 

-0.213 

0.07 

■■ 

349.821 

346.639 

-3.182 

0.91 

398.932 

392.949 

-5.983 

1.50 

B  fl 

468.320 

468.634 

0.314 

0.07 

503.548 

500.195 

-3.353 

0.67 

^B  ^B 

522.286 

521.241 

-1.045 

0.20 

560.423 

558,03 

-2.393 

0,43 

■■ 

567.995 

565.273 

-2.722 

0.48 

Table  7.  MLU  Runs  2  and  9  Frequency  Comparison  (Reproducibility) 


Figure  34.  MLU  Runs  2  and  9  Frequency  Difference  Chart  (Reproducibility) 
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MODE 

HLU21  DAMPING 

MLU9  DAMPING 

DIFFERENCE 

%  DIFF 

0.908 

2.495 

174.780 

4.054 

2.489 

38.604 

0.852 

1.074 

26.056 

0.514 

0,972 

0.458 

89.105 

5 

0.215 

0.835 

0.620 

288.372 

6 

0.587 

0.618 

0.031 

5.281 

7 

0.189 

0.674 

0.485 

256.614 

8 

0.412 

0.697 

0.285 

69.175 

9 

0.399 

0.770 

0.371 

92.982 

10 

0.227 

1.129 

0.902 

397.357 

11 

0.781 

1.354 

0.573 

73.367 

12 

0.579 

0.623 

0.044 

7.599 

13 

0.776 

0.694 

-0.082 

10.567 

14 

0.510 

0.959 

0.449 

88.039 

15 

0.299 

0.669 

0.370 

123.746 

16 

0.436 

0.767 

0.331 

75.917 

Table  9.  HLU  Run  21/MLU  Run  9  Damping  Ratio  Comparison 
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Figure  36.  HLU  Run  21/MLU  Run  9  Damping  Ratio  Difference  Chart 
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Figure  38.  MLU  Runs  2  &  9  Damping  Ratio  Difference  Chart  (Reproducibility) 
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more  to  do  with  transducer  locations  and  support  conditions  than  with  any  inherent 
difference  between  plates. 

5.  Summary  of  Results 

Even  after  accounting  for  the  reproducibility  errors,  a  measurable  difference  in 
frequencies  and  mode  shapes  is  evident  between  the  two  plates.  Because  of  the 
variability  in  damping  ratio  measurements,  the  damping  ratio  comparison  is  not 
conclusive.  If  the  frequency  and  mode  shape  results  are  assumed  to  be  a  valid 
representation  of  the  dynamic  behavior  of  the  two  plates,  it  is  clear  that  differences  in 
that  behavior  not  only  exist  but  are  quantifiable  by  reference  to  either  the  MAC  value 
or  the  frequency  shift.  These  differences  in  dynamic  behavior  can  only  be  attributed 
to  underlying  differences  in  the  stiffness,  inertial,  and  damping  properties  of  the  plates. 
In  this  regard  it  would  appear  that  the  MLU  is  either  "stiffer",  less  dense,  more 
heavily  damped,  or  some  combination  of  the  three.  (Since  the  two  plates  differ  in 
mass  by  less  than  one  gram,  the  inertial  differences  would  have  to  be  due  to  mass 
distribution.)  The  cause  of  these  differences  could  be  anything  from  less  variation  in 
fiber  volume  or  fewer  voids  in  the  MLU  plate  to  differences  in  cure  rate  between  the 
two  plates.  In  any  event,  the  significance  for  this  investigation  is  the  fact  that  the 
differences  are  identifiable. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  experimental  modal  measurement  and  finite  element  dynamic  analysis 
techniques  employed  in  this  investigation  have  been  verified  against  known  analytical 
results  for  homogeneous  isotropic  plates.  The  results  of  this  investigation  are  relevant 
to  the  following  Non-Destructive  Inspection  applications. 

1.  Quality  Assurance  and  Acceptance  Testing 

The  application  of  multi-input  multi-ouQjut  modal  measurement  and 
discrimination  as  a  Non-Destructive  Inspection  technique  for  composite  materials 
quality  assurance  applications  has  been  demonstrated  on  two  identically  configured 
composite  plates  manufactured  by  different  processes,  one  with  the  traditional  hand 
layup  method  and  one  with  machine  controlled  fiber  placement.  The  measurement 
technique  has  been  shown  to  have  both  the  reproducibility  and  resolution  to 
quantitatively  discern  the  differences  in  dynamic  responses  of  the  two  samples.  This 
quantitative  measured  difference  is  consistent  with  the  qualitative  observations  that  the 
F-18E/F  inlet  ducts  manufactured  by  machine  aided  fiber  placement  are  acoustically 
quieter  and  thus  presumably  more  heavily  damped. 

2.  Cost-Effective  Composites 

Additional  investigation  into  associating  micro-mechanical  variations  with 
stiffness  and  damping  properties  will  lead  to  refined  quantification  of  the  link  between 
internal  structural  differences  and  global  dynamics  characteristics.  Such  methodology 
may  enable  utilization  of  modal  measurement  for  defining  the  required  level  of  control 
for  materials  and  processing  in  composite  manufacturing.  For  instance,  dynamic 
response  characteristics  derived  through  modal  testing  could  be  used  as  a  quantitative 
assessment  of  the  effects  of  increasing  the  fiber  volume  variability  on  stiffness 
properties.  The  level  of  quality  control  could  then  be  set  such  that  the  desired 
stiffness  properties  were  achieved.  Such  an  approach  would  be  useful  in  the 
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manufacture  of  cost-effective  composites  for  such  non-man-safe  applications  as 
missiles  and  torpedoes. 

3.  Reliability  Assurance  and  Damage  Assessment 

The  demonstrated  reproducibility  with  removal  and  remounting  of  transducers 
opens  the  possibility  of  using  periodic  modal  measurements  (in  comparison  to  a 
benchmark  measurement  taken  when  new)  to  monitor  in-service  structural 
degradation.  This  "before  and  after"  evaluation  could  be  applied  in  both  fatigue  life 
monitoring  and  damage  assessment  situations  provided  adequate  standardization  of 
apparatus  and  procedure  could  be  achieved.  Effective  damage  assessment  would 
require  localization  of  structural  deficiencies  either  through  a  frequency  domain 
technique  such  as  that  attempted  by  Campbell  [Ref.  4],  or  by  a  comparison  of  mode 
shapes  for  each  response  coordinate.  Mode  shapes  derived  for  each  measured  spacial 
location  could  be  compared  to  those  from  the  established  benchmark  to  determine  the 
coordinates  causing  the  discrepancies  in  modal  parameters. 

4.  Field  Applications 

The  feasibility  of  adapting  modal  analysis  techniques  for  fleet  use  is  an  open 
question  that  can  only  be  resolved  with  further  refinement  of  both  apparatus  and 
procedure.  Fleet-wide  use  of  this  technique  for  NDI  would  require  standardization  of 
transducer  and  shaker  mounting  locations,  component  mounting  conditions,  and  the 
development  of  benchmark  signatures  for  each  individual  component  to  be  evaluated. 
With  the  compact  data  storage  capabilities  now  available,  compiling  a  data  bank  of 
benchmark  modal  parameters  for  critical  structural  components  would  not  appear 
problematic.  Given  that  these  requirements  can  be  achieved,  NDI  through  modal 
analysis  is  a  potentially  viable  technique  worthy  of  future  investigation. 

B.  RECOMMENDATIONS 

The  following  areas  of  investigation  are  recommended  for  future  studies  of 
NDI  using  modal  analysis. 
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1 .  Conduct  damage  localization  tests  using  before  and  after  measurements  of 
test  specimens  with  inflicted  damage.  Utilize  a  coordinate  mode  shape 
comparison  in  an  attempt  to  localize  the  damage. 

2.  Test  composite  specimens  of  various  geometries  to  verify  applicability  to 
different  component  configurations. 

3.  Use  alternate  specimen  support  conditions  to  reduce  the  variability  between 
measurements.  (A  clamped  configuration  may  be  easier  to  apply  in  the 
field.) 

4.  Mount  the  transducers  in  a  more  permanent  manner  to  further  reduce  the 
variability  between  measurements. 

5.  Refine  the  error  estimation  through  an  analysis  of  the  distribution  of  modal 
frequencies  and  damping  ratios  over  a  larger  number  of  tests.  Build  a 
probability  distribution  to  define  the  mean  and  variance  of  the  frequencies 
and  damping  ratios  for  each  specimen.  The  means  could  then  be  used  as 
more  accurate  baselines  for  comparison. 
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APPENDIX  A.  MULTIPLE  DOF  VIBRATIONS 


Consider  the  forced  2  DOF  system  modeled  by  the  dual  mass-spring-damper 
arrangement  depicted  in  Figure  A.l.  The  equations  of  motion  may  be  written  as 
follows: 

m,  X|"  +  c,  X,'  +  C2  (x/ -  X2’)  +  k,  X,  +  k2  (x,  -  X2)  =  f,  (t)  A.l 
m2  X2"  -  C2  (x,’  -  X2'  )  -  k2  (x,  -  X2  )  =  f2  (t)  A.2 

Rearranging, 

m,  X,"  +  (c,  +  C2 )  X,' -  C2  X2'  +  (k,+  k2)x,  -  k2  X2  =  f,  (t)  A.3 

m2  X2"  -  C2  X,'  +  C2  X2'  -  k2  X,  +  k2X2  =  f2  (t)  A.4 

Or,  in  matrix  form. 


Figure  A.l.  2  DOF  System 
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Now  consider  the  case  of  zero  damping  and  no  driving  force.  We  have 

[m]  {X"}  +  [k]  {x}  =  [0],  A.7 

Assuming  a  solution  x  (t)  =  A{\|/}  e‘“' ,  where  A  is  a  scalar  constant  and  {\|/}  is  a  time 
independent  amplitude  vector.  Substituting  into  equation  A.7  gives 

([k]  -  coMm])  A{\|/}  e‘“'  =  [0].  A.8 

A  non-trivial  solution  dictates  that 

det  |[k]  -  (0^  [m]  |  =  0,  A.9 

from  which  the  eigenvalues  [(0^^]  may  be  obtained.  The  eigenvectors  {\|/),  may  then 
be  found  by  substitution  of  the  eigenvalues  back  into  equation  A.8.  The  [to^^]  and  [\|/] 
matrices  represent  the  natural  frequencies  and  mode  shapes  of  the  system,  respectively. 
For  the  DEF  system  modeled  here,  we  would  have  two  natural  frequencies  and  two 
mode  shapes  corresponding  to  the  two  modes  of  vibration  (one  for  each  DOF). 

The  natural  mode  shapes  thus  described  have  the  orthogonality  property  that 
{¥}i^  [m]  {\|/}j  =0  i  #  j,  and 
(¥} M  {¥}i  =  0  i  ^  j. 

This  leads  to  the  condition 

[¥]^  [m]  [\|r]  =  [  m  J  and 
[¥]"  [k]  [¥]  =  [  k,  ]  , 

so  that  the  natural  frequency  for  each  mode  r  can  be  expressed  in  terms  of  the 
generalized  mass  and  stiffnesses  for  that  mode,  m^  and  ,  as  =  k^  /  m. 

The  orthogonality  property  allows  us  to  decouple  the  differential  equations  of 
motion  with  a  coordinate  transformation  using  the  mass-normalized  mode  shapes 
defined  as: 

[<|)f  [m]  [<!)]  =  [  I  ]  ,  A.14 

from  which  it  follows  that 

[<^]Mk]  [0]  =  [  co„M  ,  A.15 

where  [(|)]  is  the  mass-normalized  mode  shape  matrix.  The  normalized  mode  shape 
matrix  [<])]  is  then  used  as  the  transformation  matrix: 

{X}  =  m  {q},  A.16 


A.IO 

A.ll 

A.12 

A.13 
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where  {q  }  is  the  generalized  set  of  coordinates.  Substitution  into  equation  A.7 
results  in  the  following  differential  equation  of  motion: 

[m]  [<t>]  {q"}  +  [k]  [.^l  {q}  =  [0].  A.  17 

Preraultiply  by  [4>]  ^  to  obtain 

m  >][0]  {q"}  +  M  {q}  =  [0],  A.18 

or,  by  equations  A.  14  and  A.  15, 

[I]{q"}  +  [^rM{q}  =  [0].  A.19 

Equation  A.19  represents  a  set  of  uncoupled  differential  equations  that  can  be  solved 
for  {q}  and  subsequently  transformed  back  to  the  {  x  }  domain  using  equation  A.  16. 

Now  adding  a  forcing  function  {f  (t)}  =  [F]  e'"‘  to  the  right  side  of  equation 
A.  17  and  we  have 


P]  {q°}  +  {q}  =  [</»]  " 

[  F  ]  e“‘ 

A.20 

Assuming  a  solution  {q(t)}  =  [  Q]  e""*  results  in 

[F] 

A.21 

[<!)]’’  [F] 

A. 22 

To  find  the  response  in  the  frequency  domain,  we  define  the  frequency 
response  function,  or  FRF,  as  the  transfer  function  between  the  output  (displacement, 
velocity,  or  acceleration)  to  the  input  (force),  as 

Hij(a))  =  Xj/Fj  A.23 

where  X;  is  the  response  at  coordinate  i  from  the  excitation  Fj  at  coordinate  j.  For 
an  N  DOF  system,  H  (oj)  will  be  an  N  x  N  matrix  of  FRF  elements  that  constitute 
the  response  model.  Recalling  that  [Q]  =  [<^]'’  [X]  ,  equation  A.  22  gives 


_  f  <i>i  Vj 

r=l  (0)^-a)2) 


A. 24 


where  the  summation  is  made  over  all  DOF’s  (modes). 
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A.25 


For  proportional  damping  of  the  form 

[c]  =  i8[k]  +  7[m], 

where  |8  and  y  are  arbitrary  constants,  the  normal  modes  of  the  system  will  stiU 
yield  the  orthogonality  conditions  for  the  damping  matrix,  that  is, 

[c]  {\^}j  =  0  i  j,  and  A.26 

[rPY  [c]  m  =  [c,].  A.27 

Since  [c,  ]  is  also  a  diagonal  matrix,  the  system  of  differential  equations  of  motion 
remains  uncoupled  and  can  be  solved  as  an  independent  set  of  equations.  In  this  case, 
the  damped  mode  shapes  will  be  the  same  as  the  undamped  mode  shapes  and  the 
damped  natural  frequencies  are 

=  0),^  (  1  -  )  where  =  |8  oj,/  2  +  7  /  2  w, .  A.28 

The  FRF  may  be  derived  in  a  similar  manner  to  that  of  the  undamped  case  and 
is  found  to  be: 


N 


=E 


■IJ 


A. 29 


where  A'ij  =  (l>\  <!>)  ,  the  modal  constant.  An  FRF  for  a  typical  2  DOF  system  is 
shown  in  figure  A-2,  where  the  two  peaks  occur  at  the  natural  frequencies  of  the 
system  and  correspond  to  the  two  modes,  one  for  each  DOF.  FRF’s  for  most  real 
structures,  where  the  DOF’s  are  infinite,  will  exhibit  a  series  of  successive  peaks, 
some  of  which  are  more  dominant  than  others,  but  each  of  which  correspond  to  a 
particular  vibrational  mode  of  the  structure.  The  damping  at  a  particular  mode  is 
seen  as  a  blunting  of  the  peaks  in  the  FRF.  [Ref.  5,  7,  14] 
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APPENDIX  B.  MODAL  PARAMETER  EXTRACTION 


The  following  describes  the  process  used  to  extract  the  modal  parameters  from 
the  set  of  measured  FRF's. 

The  first  step  of  the  extraction  procedure  was  to  initiate  a  modal  parameter 
table  within  IDEAS  Modal  Analysis  by  determining  the  resonant  modes  and  an 
estimate  of  the  frequencies.  A  Multivariate  Mode  Indicator  Function  (MMIF)  routine 
was  used  to  identify  the  modes  and  their  frequencies.  The  MMIF  minimizes  the  ratio 
of  the  kinetic  energies  of  the  in-phase  response  (real  part)  to  the  total  response  for  the 
entire  set  of  FRF's  to  identify  the  normal  modes.  This  idea  is  based  on  the  fact  that  at 
the  normal  modes,  the  imaginary  of  the  response  reaches  a  maximum  and  the  real  part 
tends  to  zero.  The  Mode  Indicator  Function  (MIF)  for  a  single  FRF  is  characterized 
by  the  function: 


E  W 


B.l 


where  H  is  the  FRF.  The  MMIF  is  an  extension  of  this  concept,  but  it  finds  the  force 
vector  that  produces  the  smallest  ratio  of  real  response  to  total  response. 


F'^AF  ^  ^ 
F'^{A+B)F 


B.2 


where 

FT  A  F  =  F^  (H^„  M  F 

F"  B  F  =  F^  M  )  F 

a  =  a  diagonal  matrix  of  the  minima. 

The  minima  are  found  by  solving  the  eigenvalue  problem 

AF=  (A  +  B)Fa,  B.3 

and  the  eigenvalues  can  then  be  plotted  as  a  function  of  frequency  to  reveal  the  real 
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normal  modes.  Localized  or  lower  amplitude  modes  are  revealed  by  the  MMIF  since 
the  maximum  amplitude  is  normalized  to  one. 

The  MMIF  has  a  maximum  value  of  one  with  downward  spikes  at  the 
frecjuencies  corresponding  to  the  normal  modes.  Figure  B.l  shows  the  MMIF  for  the 
MLU  test  FRF's.  A  "search  valleys"  function  picks  the  modal  frequencies  from  the 
MMIF  with  a  user  defined  threshold  for  maximum  values  of  the  spikes.  The  end 
product  of  this  step  was  a  parameter  table  listing  the  frequencies  of  the  identified 
modes  for  the  set  of  measured  FRF's. 

The  Polyreference  complex  exponential  curve  fitting  method  was  then  used  to 
estimate  the  damping  ratios  and  refine  the  resonant  frequency  estimates.  The 
Polyreference  technique  is  a  time  domain  algorithm  that  simultaneously  processes 
multiple  FRF's  to  obtain  a  global  least  squares  estimate  of  the  modal  parameters.  The 
method  was  developed  by  Void  et.  al.  [Ref  7].  Essentially,  it  fits  decaying 
exponentials  to  impulse  response  functions  through  an  inverse  Fourier  transform  of  the 
FRF's.  A  summary  of  the  theory  behind  the  technique  follows. 

The  response  to  a  unit  impulse  at  an  exciter  location  can  be  obtained  from  the 
inverse  Fourier  transform  of  the  appropriate  FRF.  In  matrix  form,  it  is: 

{X.  (t,  )}  =  [W][  ]  {A,j  }  B.4 

where 

Xj  (tfc )  =  an  1  dimensional  vector  representing  the  response  at  i  due  to  a  force  at  j, 
where  1  is  the  total  number  of  excitation  locations 
W  =  an  1  X  2N  matrix  representing  the  ratios  of  the  eigenvectors  of  the  different 
exciter  locations 

U'‘  =  a  2N  X  2N  diagonal  matrix  of  diagonal  elements  e*r'‘^ 

Ajj  =  a  2N  dimensional  vector  of  the  residues  for  response  i  and  exciter  j 
N  =  number  of  modes 
T  =  time  subinterval 
k  =  1,2,3, 4,.. .2N 
A  matrix  G  may  be  defined  as 
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Figure  B.l.  Multivariate  Mode  Indicator  Function  for  MLU  Test  Run  2 
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where  p  is  the  smallest  integer  such  that  p  1  >  2N. 

Then  there  exists  a  full  rank  matrix  Q  such  that 

[Q]  [G]  =  0. 

Q  can  be  partitioned  into  square  submatrices 

[Q]  =  [  [B(p)]  [B(p-l)]  . [B(l)]  [B(0)]  ] 

Equation  B.6  can  then  be  rewritten  as 


'f^B{s)WUP'^  =  0  B.8 

S=0 

The  coefficients  B  are  found  by  multiplying  each  coefficient  by  the  response  given  by 
equation  B.4  for  the  specific  time  t^  resulting  in 


This  equation  is  solved  using  a  least  squares  method  to  obtain  the  coefficients  B.  A 
byproduct  of  this  procedure  is  the  lower  triangular  correlation  matrix  containing  the 
residual  error  terms  used  to  determine  the  number  of  poles  required  for  the  analysis. 
An  error  chart  and  pole  stability  diagram  are  displayed  when  the  Polyreference 
technique  is  implemented  in  IDEAS,  allowing  the  selection  of  the  appropriate  number 
of  poles  when  building  the  parameter  table. 

To  determine  the  resonance  frequencies  and  damping  values,  the  method  finds 
the  zeroes  of  the  matrix  polynomial  enclosed  in  parenthesis  in  the  following 
expression 

Equation  B.IO  was  derived  by  multiplying  equation  B.8  above  by  e, ,  the  vector  of 
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B.  10 


P 

(J^B(s)  =  0 


zeroes  except  for  a  one  in  position  r.  Once  [U]  is  known,  the  frequencies  and 
damping  values  follow  from 

U  ==  ^  j?»dr)T  B  1  ] 

and  ln|U,|/T=  -Cr  ,  Z  U,  /  T  =  w*  =  (O,  (1  -  Cr  )'''•  B.12 

The  eigenvectors  may  then  be  found  in  the  null  space  associated  with  the  zeroes 
of  the  matrix  polynomial. 

The  Frequency  Polyreference  technique  was  used  to  generate  the  residues  and 
mode  shapes  after  the  natural  frequencies  and  damping  ratios  had  been  determined. 
The  Frequency  Polyreference  method  has  the  advantage  of  operating  in  the  frequency 
domain,  permitting  the  use  of  the  residual  stiffness  and  inertia  terms  from  modes 
outside  the  frequency  range  of  interest  in  calculating  the  mode  shapes.  The  relation 
between  the  FRP's  and  the  residues  may  be  written  in  matrix  form  as  follows: 


'Aij 

= 

w'u^ 

WU^ 

Ail 

where 

Hj  (IcOl  )  =  an  1  dimensional  vector  consisting  of  the  value  of  the  FRP  at  the 
frequency  (Ol  for  a  response  at  i  and  an  exciter  at  j 
W  =  the  same  matrix  as  above  except  that  it  has  been  scaled  so  that  the  row 
corresponding  to  the  exciter  has  is  set  to  one 
U„  =  a  2N  X  2N  diagonal  matrix  with  the  r'**  diagonal  element  equal  to  1  /  (s  -  s„ ) 
L,  H  =  frequency  range  of  interest. 

Equation  B.14  can  be  written  as 
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[H]  =  [P][A],  B.14 

and  the  least  squares  solution  is  found  by  the  normal  equations  method,  as  follows: 

[Pf  [H]  =  [[Pf  [P]][A]  B.15 

[[Pf  [P]]‘  [[Pf  [H]]  =  [A].  B.16 

The  residual  compliance  and  inertance  terms  can  be  added  into  the  [  W  ]  and  [  U  ] 
matrices  prior  to  solving  the  equation  for  the  residues. 

With  the  natural  frequencies,  damping  ratios,  and  residues  calculated,  the 
parameter  table  was  complete.  The  residues  are  displayed  as  an  amplitude  value  and 
a  phase  angle.  The  analytic  function  representing  these  parameters  can  be  generated 
and  compared  to  the  measured  FRF’s  to  qualitatively  check  the  fit;  two  representative 
samples  were  shown  in  Figures  17  and  18.  If  the  fit  is  not  satisfactory,  the  process 
can  be  repeated  using  a  different  number  of  poles  or  with  the  lower  valued  MCF 
modes  filtered  out. 

The  final  step  in  the  process  is  that  of  generating  the  mode  shapes.  This  is 
accomplished  by  the  Polyreference  algorithm  from  the  residues,  damping  ratios,  and 
resonant  frequencies,  since 

A'  =  A^„  +  A^  B.17 


A  ^ 

“re  r 

4r<*>r 


B.18 


'  i/n  r  2 

Cr  Wr 


B.19 


In  the  case  of  real  normal  modes,  A^„  =  0  so  that  =  0  and 


B.20 


The  mode  shapes  are  stored  as  a  table  of  coefficients  representing  the  deformation  at 
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each  node.  Real  or  complex  mode  shapes  can  be  generated  using  either  3  or  6  DOF's. 
The  mode  shapes  from  two  different  tests  or  from  a  finite  element  model  solution  and 
a  test  can  then  be  compared  using  Modal  Assurance  Criteria  or  by  graphical  portrayal 
of  the  deformations.  [Ref.  5,  7,  15] 
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APPENDIX  C.  MODAL  ASSURANCE  CRITERIA 


The  following  presents  the  background  behind  the  Modal  Assurance  Criteria 
(MAC)  correlation  technique. 

The  MAC  provides  a  means  to  compare  two  mode  shapes,  whether  they  were 
derived  from  two  separate  tests,  from  a  test  and  a  Finite  Element  solution,  or  from 
the  same  test  using  different  extraction  techniques.  MAC  is  a  scalar  value  between 
zero  and  one  that  represents  the  degree  of  correlation  between  the  two  mode  shapes. 
A  value  of  one  indicates  good  correlation  where  the  mode  shapes  differ  only  by  a 
scale  factor;  a  value  of  zero  indicates  poor  correlation. 

The  MAC  is  a  least  squares  estimation  technique  defined  by  the  following 
expression: 


MAC  = 


^»pi(r)cpi(r)  52  <p2(r)<P2(r) 


C.l 


where  i/'i  (  r  )  =  mode  shape  for  rth  degree  of  freedom 
^j’  (  r  )  =  mode  shape  complex  conjugate 
m  =  number  of  degrees  of  freedom 

The  MAC  must  be  used  with  some  caution,  since  there  are  other  factors  which 
can  influence  its  value.  For  instance,  the  MAC  may  assume  a  value  near  zero  for  the 
following  reasons: 

1.  The  tested  system  is  nonstationary  due  to  changes  in  mass,  stiffness,  or 
damping  properties  during  the  test. 

2.  The  system  is  non-linear. 

3.  The  reference  mode  shape  contains  noise. 

4.  The  method  used  for  parameter  extraction  is  invalid  for  the  measured  data. 

5.  The  mode  shapes  are  linearly  independent. 
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Likewise,  the  MAC  can  take  on  a  value  near  unity  for  the  following  reasons: 

1.  The  number  of  responses  is  not  sufficient  to  distinguish  between  mode 
shapes. 

2.  Unmeasured  forces  were  responsible  for  the  mode  shapes. 

3.  The  mode  shapes  are  mainly  coherent  noise, 

4.  The  mode  shapes  represent  the  same  motion  differing  only  by  a  scale 
factor. 

Obviously,  it  is  the  last  reason  in  each  of  these  lists  that  is  desired;  if  the  others  can 
be  eliminated  the  MAC  can  be  inteipreted  as  an  approximate  orthogonality 
calculation.  [Ref.  5,7] 
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APPENDIX  D.  TEST  SPECIMEN  MATERIAL  PROPERTIES 


The  laminates  used  in  the  experiment  were  constructed  from  Hercules  IM7 
carbon  fiber  and  8552  epoxy  matrix.  IM-7  is  a  continuous,  intermediate-modulus, 
polyacrylonitrile  (PAN)  based  fiber  that  combines  high  tensile  strength  and  modulus 
with  good  shear  strength.  The  laminate  had  16  plys,  each  approximately  .0055  inches 
thick,  in  a  symmetrical  -1-45  /  0  /  -45  /  -b90  layup  sequence.  Listed  below  are  the 
specific  mechanical  properties  of  the  ply  and  the  laminate  [Ref.  16]. 


Property 

Ply 

Laminate 

Axial  Tensile  Modulus,  Ex  ,  Msi 

25 

9  4* 

Transverse  Tensile  Modulus,  Ey  ,  Msi 

1.39 

9.4* 

Axial  Tensile  Strength,  X  ,  ksi 

400 

— 

Axial  Compression  Strength,  X’,  ksi 

240 

— 

Shear  Strength,  S,  ksi 

14.5 

— 

Density, p  ,  Ib/in^ 

.06 

.056t 

Shear  Modulus,  Gxy  ,  Msi 

.71 

3.6* 

Shear  Modulus,  Gyz  ,  Msi 

.07* 

.13* 

Poisson  Ratio,  vxy 

.31 

.32* 

Poisson  Ratio,  vyz 

.34* 

.32* 

Fiber  Volume,  Vf 

.62 

.62 

t-Measured  ^-Estimated  from  micromechanics  relations. 

Table  D.  1  Material  Properties  of  Test  Specimens 
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Figure  D.2  Machine  Layup  Plate  (Front  and  Back) 
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APPENDK  E.  MECHANICAL  BEHAVIOR  OF  LAMINATES 


The  mechanical  behavior  of  a  laminate  is  governed  by  the  constitutive  matrix 
Cij  ,  which  relates  the  applied  stress  to  strains  and  curvatures.  The  constitutive 
matrix  is  made  up  of  a  number  of  submatrices,  as  seen  in  the  following  expressions 
that  relate  the  forces  and  moments  to  the  displacements  [Ref.  10,  p.  154]. 


f^xl 

All 

Ai2 

VO 

H 

fill 

S12 

S16 

Ny\ 

Ai2 

A22 

A26 

1  ^0  1 
\  y  \ 

►  + 

S12 

B22 

S26 

K3,  . 

KJ 

Ai6 

A26 

^66. 

[yvI 

S16 

^26 

^66. 

f^xl 

’^11 

Bi, 

’^11 

^12 

\^x' 

My\ 

►  = 

B22 

B2e 

>  + 

D^2 

B22 

^26 

- 

K\ 

.^16 

S26 

^66. 

[yvJ 

.^16 

^26 

^66. 

where 

N,  M  =  resultant  forces  and  moments 
e  =  strain 
7  =  shear  strain 
K  =  curvature 

A  =  extensional  stiffness  matrix 
B  =  coupling  stiffness  matrix 
D  =  bending  stiffness  matrix. 

For  regular  symmetric  angle-ply  laminates,  the  plies  are  of  equal  thickness 
and  are  arranged  such  that  the  orientation  angles  of  adjacent  plies  are  opposite  in  sign 
and  symmetrically  layered  about  the  middle  of  the  laminate.  The  test  specimens  used 
in  this  investigation  were  a  special  form  of  this  type  of  laminate,  since  the  opposing 
angle  plies  of  +/-  45  degrees  were  separated  by  layers  oriented  at  0  or  90  degrees. 

For  symmetric  laminates  the  terms  of  the  coupling  matrix,  B,  are  zero, 
eliminating  any  coupling  between  bending  and  extension.  However,  the  presence  of 
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the  terms  Dj^  and  Djg  introduces  twist-curvature  coupling  for  these  types  of 
laminates.  This  t3Tje  of  coupling  has  been  shown  to  reduce  the  effective  stiffness  of 
the  laminate  with  a  resulting  increase  in  deflections  and  a  decrease  in  vibration 
frequencies.  The  effect  of  this  coupling  on  deflections  and  vibration  frequencies  dies 
out  as  the  number  of  ply  layers  is  increased. 
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APPENDIX  F.  SIMULTANEOUS  VECTOR  ITERATION 


The  following  describes  the  method  of  Simultaneous  Vector  Iteration  (  SVI ) 
as  applied  to  structural  dynamics  problems  [  Ref.  6  and  13]. 

SVI  is  two  step  inverse  power  iteration  method  that  solves  the  eigenvalue 
problem 

[k]{<f>}  =  X[m]{<A}  F.l 

where  X  is  the  square  of  the  natural  frequencies  and  {<^>}  is  the  mode  shape  vector. 
Each  iteration  in  the  algorithm  consists  of  two  steps; 

1 .  Inverse  power: 

]  =  [  k  ]-'  [  m  ]  [<^  ]p  F.2 

2.  Subspace  orthogonalization: 

[<t>  ]p.,  =  ]  [Q]  F.3 

where 

[ypf  [k][^][Q]  =  [m]  [t^][Q][X]  F.4 

and 

[<!>]  =  mode  shape  matrix 
[  m  ],  [  k  ]  =  mass  and  stiffness  matrices 
[  Q  ]>  [  ^  ]  =  eigenvectors  and  eigenvalues  of  equation  F.4 
p  =  iteration  counter 

As  the  iteration  proceeds,  [  ^  ]  and  [  X  ]  converge  to  the  mode  shapes  and  the  square 
of  the  natural  frequencies  of  the  system. 
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APPENDIX  G.  FEE’S  THROUGH  RANDOM  EXCITATION 


For  a  stationary  ergodic  random  process  f  ( t ),  the  autocorrelation  function  is 
defined  as: 


T/2 

Rf^it)  =  E{f(t)  £(t+T)}  =  ^  jf  fit)  f{t+x)  dx  G.l 

-T/2 

where  T  is  the  total  time  of  the  sample  and  E{  }  denotes  the  expected  value  of  the 
product  of  f  at  time  t  and  t  +  t.  The  autocorrelation  function,  unlike  the  random 
process  it  was  derived  from,  meets  the  Dirichlet  criteria  for  discontinuities  and 
therefore  the  Fourier  transform  may  by  obtained  as  follows: 


~00 


G.2 


where  Sff  is  the  auto-spectral  density  function.  The  auto-spectral  density  function 
provides  frequency  domain  characterization  of  the  random  process. 

Similarly,  for  a  pair  of  functions  f(t)  and  x(t),  the  cross  correlation  and  cross 
spectral  density  are  defined  as: 


T/2 

=  £Gc(t)  f(t+T)}  =  1  J  xit)  fit+x)  dx  G.3 

-T/2 

and 


Now  consider  an  input  process  f(t)  and  a  response  x(t).  In  order  to  find  a 
relation  between  the  input  and  response  spectral  densities  we  introduce  the  impulse 
response  function  g(t)  and  the  dummy  variables  X,  and  X2  into  the  convolution 
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integrals  as  follows: 


x{t)  =  j  f(t-X^)  dX^ 


x{t+x)  =  Jsr(A,2)  r  ( t+x-Xj)  c?A,2 

— oo 

Using  equation  G.l  above  and  the  fact  that  E{g(x)}  =  g(x),  we  can  write 

«  oo 

R^ix)  =E(f  sr(X^)  f(t-X^)  dX,  f  ff(X2)f(t^x-X,)dX,}  G.7 


=  JJ  g[X^)  g{X^)  E{f(t-X^)  /( t+T-Xj )}  dXidAj 
—  00  —  00 

An  ergodic  process  is  stationary;  hence, 

E{f{t-X^)  f{t+x-X2))=E{f(t)  f{t+x+X^-X2)) 


=  Rff{x+X^-X2) 


G.IO 


which  is  the  autocorrelation  of  the  excitation  process.  Substituting  into  equation  G.8 


gives 


R^x ^x)  ~  J'  J  Xj^)  g ( X2 )  ( t+T  X2 ) )  dX-^dX2 


G.ll 


Substituting  this  equation  into  the  expression  for  the  autospectra  of  the  response  yields 

00  00  OQ  00 

e-^^^dx  =  J  [|  J  g(X^)  g(X2)  Rff(t+x-X2)  dX^dX2 


G.12 
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Since 


G.13 


upon  substitution  into  equation  G.12  we  obtain 


5^(0))  =  ga^)ga^)  [ f  5^,^ (a>)  do)]  dX.dX^ ) 


G.14 


jf  J5^f(a))  [jf  gr(Xi)e^"^dXi  |  sr(^2)  dT 


G.15 


C»  OD 


G.16 


OO  00 

=  J  e"^“^  Js'f:/ (i))  |fj(  w)  P  e^“^d(«)}  dT  G.17 

—00  -CO 

where  H(-a))  is  the  complex  conjugate  of  the  frequency  response  H(&)).  Taking  the 
inverse  Fourier  transform  of  equation  G.17  we  have 


t'S'xx]  e-^dw=^/5ff(a))  \H(<o)  |2e^"^d(0 


G.18 


5^(0))  =  |H(<i>)  p  (w) 


G.19 


This  expression  relates  the  spectral  densities  of  the  input  and  response  through  the 
transfer  function  H(a)).  In  a  similar  fashion  the  expressions  for  the  cross  spectra  are 
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s^{<x>)  =  H((ji)  S^f((ji) 


G.20 


Sf^{(o)  =Hi<o)  5ff(o>)  G.21 

Equations  G.19,  G.20  and  G.21  make  it  possible  to  obtain  the  FRF  from  a 
random  excitation  in  at  least  three  different  manners.  This  provides  a  convenient 
check  on  the  data,  since  each  method  should  give  the  same  FRF.  The  ratio  of  the 
various  transfer  functions  is  known  as  the  coherence,  which  ideally  should  be  equal  to 
1.0  if  only  one  excitation  force  is  used. 

In  the  IDEAS  software  and  elsewhere  the  standard  nomenclature  for  the 
transfer  functions  is  Hj  and  Hj  for  the  functions  H  in  equations  G.20  and  G.21, 
respectively.  IDEAS  has  the  capability  to  produce  the  FRF  through  both  of  these 
methods.  In  addition,  IDEAS  can  also  calculate  the  FRF  using  the  optimal  scaling 
method,  which  formulates  a  noise  spectral  matrix  to  account  for  truncation  errors. 
[Ref.  5,  7,  14] 
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APPENDIX  H.  DATA  ACQUISITION  SETTINGS 


Sampling  Conditions 


Spectral  Lines  :  801  Delta  f  :  1.25  Hz 

Maximum  Freq  ;  1000.0  Hz 

Minimum  Freq  ;  0.0  Hz 

Frame  size  :  2048  Frame  Length  :  0.8  sec 

Delta  t  :  3.91E-04  sec  Sampling  Freq  :  2560.0  Hz 

Replay  Ratio  :  1.0 

Trigger  Conditions 


Trigger  Method  :  Free  Run  Delay  Type  ;  Pre  Trigger 

Trigger  Source  ;  Channel  Delay  Samples  :  204 

Input  Channel  :  1  Delay  Time  :  0.07969 

Level  Percent  :  5.0  Delay  Percent  :  10.0 

Slope  :  Positive 

Trigger  Bell  :  On 

Autorange  Conditions 


Autorange  Before  Preview  :  On  Percent  Frame  :  100.0 

Autorange  Before  Acquire  :  On  Percent  Overhead  :  30.0 

Window  Conditions 


Ref /Res  Window  :  Hanning  Broad 

Averaging  Conditions 


Average  Method  :  Stable 

Exponential  Average  Constant  :  5 

Frames  per  Average  :  20 

Overlap  Percent  :  0.0 

Frame  Acceptance  :  Automatic 

Measurement  Conditions 


Acquisition  Results 

:  FRF 

Auto  Switch 

;  Off 

Clear  Lower 

:  0.0 

Clear  Upper 

:  1000.0 

Coherence  Switch 

:  Off 

Normalization  ; 

Units"2/Hz 

Ref  Coherence  Switch 

:  On 

Amplitude  Units 

:  RMS 

FRF  Switch 

;  On 

FRF  Method  : 

Hl-Gyx/Gxx 

Test  Log  Switch 

:  On 

Log  Prefix 

:  TEST9 

Function  Logging  Switch 

;  Off 

Description 

:  MLU 

PLATE  TEST  9  /  2C 

1  MAR  95 

Real  Time  Display  Conditions 


Number  of  Display  Channels 
Display  Units 
Background  Grid 
Range  Indicators 
Upper  Limit 

Preview  Sampling  Factor 
Acguisition  Monitor 


8 

Engineering  Units  (EU) 

None 

On 

90.0  Lower  Limit  : 
1.0 

Time  and  Windowed 


Waterfall  Options 


Function  Height  Percent  :  100.0 
Number  of  Functions  :  30 

Start  Amplitude  Percent  :  100.0 
Hidden  Line  switch  :  On 


Overall  Shutdown  Conditions 


Overall  Shutdown  Switch  :  Off 


Signal  Sources 


Signal  Sources 
Signal  Type 
Duration  (%) 

Off  Time  (%) 
Lower  Freq  Bound 
Upper  Freq  Bound 

Source 


:  Off 

:  Burst  Random 
:  30.0 
:  70.0 
:  0.0  Hz 
:  1000.0  Hz 

Level 


1- Sl  0.09899 

2- S2  0.09899 


20.0 
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APPENDIX  I.  CALIBRATION  PROCEDURES 


A.  SIGNAL  ANALYZER  CALIBRATION 

The  Zonic  System  7000  signal  analyser  was  calibrated  automatically  upon 
initialization  from  IDEAS  Test  and  whenever  signal  processing  settings  were  changed. 
The  System  7000  calibrates  itself  through  an  internally  generated  signal  that  checks 
gain  settings,  filtering  and  phase. 


B.  TRANSDUCER  CALIBRATION 


Each  transducer  was  individually  calibrated  using  a  mass  calibration  procedure 
whereby  a  block  of  aluminum  of  known  mass  was  connected  to  one  of  the  shakers 
through  a  force  transducer.  Each  accelerometer  was  then  placed  opposite  the  force 
transducer  one  at  a  time  and  a  random  excitation  was  applied.  The  result  was  a  flat 
response  at  a  level  corresponding  to  the  inverse  of  the  mass,  since 

F  =  ma  I.l 

a  /  F  =  1  /  m  1.2 

and  the  FRF  we  measured  is  simply  a  /  F.  A  scale  factor  was  applied  to  each 
transducer  channel  to  adjust  the  input  voltage  level  to  meet  the  below  criteria: 


1.3 


For  example,  a  sample  calibration  FRF  before  application  of  the  scale  factor  is 
shown  in  Figure  I.l.  The  flat  portion  of  the  curve  is  at  a  level  corresponding  to 
approximately  1791  in  /  Ibf  -sec^  .  The  mass  of  the  calibration  block  and  both 
transducers  was.0006115  Ibf  -sec^  /  in.  ,  so  that 


_  1  _  1  *  1 
m  Vg  .0006115  1791 


0.9131 


1.4 


The  scale  factor  so  obtained  was  then  applied  to  the  channel  in  question  through  the 
channel  set  up  feature  in  IDEAS  Test.  Figure  1.2  shows  the  calibration  FRF  after 
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Acce I  era  1 1  on  Phase 


application  of  the  scale  factor,  where  the  level  is  within  one  percent  of  the  l/m  value 
of  1635.3  in/lbj-sec^2.  This  procedure  was  repeated  for  each  transducer. 

C.  SHAKER  CHECK  OUT 

Prior  to  the  initial  test,  each  shaker  was  checked  using  the  procedure  detailed  in 
the  Ling  Electronics  Model  207  Operation  Manual.  The  procedure  involves 
connecting  an  oscilloscope  to  the  output  of  an  accelerometer  attached  to  the  shaker's 
mounting  surface  and  running  a  swept  sine  signal  through  the  shaker.  The 
accelerometer  output  is  then  monitored  for  any  signal  anomalies.  The  procedure  was 
conducted  from  0  to  2000  Hz  with  no  discrepancies  noted  for  either  shaker. 

D.  TRANSDUCER  SUITABILITY  CHECK 

To  ensure  the  available  transducers  were  suitable  for  the  selected  frequency 
range,  the  low  frequency  response  was  calculated  from  the  Discharge  Time  Constant 
(DC)  according  to  the  PCB  accelerometer  manual.  Specifically,  the  3  DB  down 
frequency  is  approximately 

f,  =  .16  /  DC  1.5 

for  the  Model  336C34  transducers.  The  DC  of  our  transducers  was  50  sec,  or  the 
low  frequency  cut  off  was  .0032  Hz,  well  below  the  first  bending  mode  of  our 
structure. 

The  upper  frequency  limitation  was  governed  by  the  resonance  frequency  of 
the  transducers,  which  was  listed  as  70  KHz.  The  PCB  manual  states  that  the 
transducers  should  only  be  used  for  applications  within  the  lower  20  percent  of  the 
resonance  frequency,  or  below  14  KHz.  This  frequency  was  well  above  our  range  of 
interest. 
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APPENDIX  J.  SAMPLE  DATA  OUTPUT 


MODE 

FREQUENCY 

(HERTZ) 

DAMPING 

(%) 

AMPLITUDE 

PHASE 

(RAD) 

MCF 

1 

64.084 

0.908 

8.2953E+04 

1.571 

0.994 

2 

82.878 

4.054 

1.6422E+05 

1.571 

0.982 

3 

110.628 

0.852 

3.4760E+04 

-1.571 

0.995 

4 

133.179 

0.514 

2.1945E+05 

-1.571 

0.999 

5 

159.599 

0.215 

2.5971E+05 

1.571 

0.998 

6 

221.542 

0.587 

1.3544E+05 

1.571 

0.990 

7 

238.048 

0.189 

5.2726E+05 

-1.571 

0.999 

8 

288.421 

0.412 

7.3539E+04 

1.571 

0.994 

9 

306.315 

0.399 

1.2519E+05 

-1.571 

0.994 

10 

343.383 

0.227 

4.7727E+05 

1.571 

0.997 

11 

390.374 

0.781 

7.7894E+05 

1.571 

0.984 

12 

465.919 

0.579 

1.7857E+05 

-1.571 

0.994 

13 

494.100 

0.776 

1.3248E+05 

1.571 

0.997 

14 

505.978 

0.510 

1.8766E+05 

-1.571 

0.999 

15 

553.979 

0.299 

1.6830E+06 

1.571 

1.000 

16 

562.869 

0.436 

5.9721E+05 

1.571 

0.999 

Table  J.l.  Modal  Parameters  Derived  by  Polyreference 
Technique,  HLU  Test 


MODE 

FREQUENCY 

(HERTZ) 

DAMPING 

(%) 

AMPLITUDE 

PHASE 

(RAD) 

MCF 

1 

67.955 

3.081 

5.8694E+04 

-1.571 

0.998 

2 

86.887 

4.407 

2.1013E+05 

-1.571 

0.987 

3 

114.458 

1.621 

1.0674E+05 

1.571 

1.000 

4 

138.654 

1.591 

1.1914E+05 

1.571 

0.997 

5 

163.659 

1.318 

3.7298E+05 

1.571 

0.999 

6 

221.525 

0.746 

6.1483E+05 

1.571 

0.999 

7 

242.358 

0.902 

5.5301E+05 

1.571 

1.000 

8 

291.394 

0.810 

1.6802E+05 

1.571 

0.999 

9 

304.421 

0.925 

1.3300E+05 

1.571 

0.998 

10 

346.639 

0.975 

6.5222E+05 

-1.571 

0.997 

11 

392.949 

2.480 

7.0102E+05 

1.571 

0.958 

12 

468.634 

0.905 

3.8168E+05 

1.571 

0.998 

13 

500.195 

0.906 

3.0744E+05 

1.571 

0.977 

14 

521.241 

1.098 

5.6214E+05 

1.571 

1.000 

15 

558.030 

0.621 

8.0497E+05 

1.571 

0.999 

16 

565.273 

0.920 

1.4245E+06 

1.571 

0.999 

Table  J.2.  Modal  Parameters  Derived  by  Polyreference 
Technique ,  MLU  Test 
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Mode 

HLU  Freq 

MLU  Freq 

%  Change 

MAC 

FE  Freq 

1 

63.914 

67.955 

6.323 

0.981 

63.12 

2 

83.274 

86.887 

4.339 

0.960 

88.50 

3 

110.916 

114.458 

3.193 

0.942 

112.27 

4 

133.177 

138.654 

4.113 

0.987 

148.74 

5 

159.609 

163.659 

2.538 

0.978 

157.16 

6 

221.529 

221.525 

-0.002 

0.966 

247.84 

7 

238.054 

242.358 

1.808 

0.992 

251.21 

8 

288.421 

291.394 

1.031 

0.960 

300.40 

9 

306.314 

304.421 

-0.618 

0.804 

321.46 

10 

343.377 

346.639 

0.950 

0.954 

347.88 

11 

390.334 

392.949 

0.670 

0.789 

403.10 

12 

465.959 

468.634 

0.574 

0.664 

13 

494.085 

500.195 

1.237 

0.056 

493.75 

14 

505.957 

521.241 

3.021 

0.520 

526.34 

15 

553.995 

558.030 

0.728 

0.840 

551.64 

16 

562.851 

565.273 

0.430 

0.723 

563.87 

Table  J.3.  Modal  Assurance  Criteria 
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